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are presently available for the IBM—PC and IBM mainframe. The Apple Macintosh 
offers just a one-parameter root locus method capability. MacRootLocus supports 
both one— and two—parameter root locus methods on the Apple Macintosh. It is 
written in the computer language Turbo Pascal which is the native language of the 
Apple Macintosh and is designed with the same user—friendliness and standard 
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I. MACROOTLOCUS SYSTEM 


A. INTRODUCTION 

MacRootLocus system is a program written for the Apple Macintosh 
computer. It was designed to allow a user to analyze and design linear feedback 
control systems. It is written in the computer language Turbo Pascal which is the 
native language of the Apple Macintosh and designed with the same 
user—friendliness and standard interface philosophy the Macintosh was designed for. 

The Macintosh gets most of its input from the user through the 'mouse’ which 
is a small cigarette pack—sized control that is moved across the table much like a 
pencil across paper to move the cursor or pointer on the computer screen. Rather 
than typing in commands like the IBM, the Macintosh lets you select the function 
you want performed with the mouse. Since the commands are not typed in, the 
commands do not have to be remembered as with other computers. It is the 
user—friendliness that sets MacRootLocus apart from other system analysis 
programs. Prior computer experience is not required to use MacRootLocus. A few 
minutes to learn how to use the mouse and pull down menus is all that 1s needed. 
All dialogs are easy to use and there are on—line Help menus. So the first-time user 
can get desired results without using trial and error. 

Apple Macintosh models supported by MacRootLocus are the Macintosh SE 
and the Macintosh II system. 

MacRootLocus was tested with several examples and is now available to any 
user on the Naval Postgraduate School Controls Laboratory computers under the 


icon of MacRoot Locus system. 


B. MATHEMATICAL CONCEPT 
The root locus method is a graphical technique for determining the roots of the 
closed—loop characteristic equation of a system as a function of the static gain. 


Consider the general feedback control system, as shown in Figure 1.1. 





Figure 1.1 Closed—loop for Control System 


In order to find the roots of the characteristic equation, it is required that 


1 + G(s)H(s) = 0. (1.1) 


Of course, Eq. (1.1) may be rewritten as 


G(s)H(s) = =I. (152) 


since s is a complex variable, Eq. (1.2) may be rewritten in parametric form 


|G(s)H(s)| = 1 (le) 





L G(s)H(s) = (2k-1)z (1.4) 


where k is an integer. 

For a specific value of s to be a root of the characteristic equation, it must 
satisfy both Eqs. (1.3) and (1.4). Since the roots are those values of s that satisfy 
both equations, then the root points are points where these curves intersect. 

The original development of the root locus method was concerned with the 
determination of the locus of roots of the characteristic equation as the system gain, 
K, is varied from zero to infinity. It appears that the root locus method is a single 
parameter method; fortunately it can be readily extended to the investigation of two 
or more parameters. 


The characteristic equation of a dynamic system may be written as 


sea otter, aS a So. OES) 
n n-1 1 0 
Clearly, the effect of the coefficient a; may be ascertained from the root locus 


equation 


A.D 
SL SS (1.6) 
Pe eee Gu) ee he acy GCE. a 
n n-] 2 0 


The parameter of interest, A (in MacRootLous), can be isolated as 


a S"+a Pr t....¢(a —A)S™O4 ASDA + |, 
n n- n-q 


+tasSt+aos=0. (27) 
1 0 


For example, a third—order equation of interest might be 
S? + (3 + A)S? + 38+6=0. (1.8) 


In order to ascertain the effect of the parameter A, we isolate the parameter and 


rewrite the equation in root locus form as shown in the following steps: 
S3 + 382 + AS? + 35 +6=0 (1.9) 


2 
en (1.10) 


Go aS Caen 


Then, to determine the effect of two parameters, we must repeat the root locus 
approach twice. Thus for a characteristic equation with two variable parameters, A 


and B (in Mac RootLocus), 


a S™-a So) la — Ao ee ee 
n ne} n-q 


+(a —B)S"'+4+ BST? 4 ..4 aSta =0. (115 


Nek 


The two variable parameters have been isolated and the effect of A will be 
determined, followed by the determination of the effect of B. For example, for a 


certain third—order characteristic equation with A and B as parameters, we obtain 


SoS) Boe (1 


In this particular case, the parameters appear as the coefficients of the characteristic 
equation. The effect of varying B from zero to infinity is determined from the root 
locus equation 


| 3. ee ee (1.13) 


S°4+8? + A 
One notes that the denominator of Eq (1.13) is the characteristic equation of the 


system with B = 0. Therefore, one first evaluates the effect of varying A from zero 


to infinity by utilizing the equation 


S°+S7+A=0 ile 


rewritten as 


[eens Se ag (1.15) 


S7(S + 1) 


where B has been set equal to zero in Eq. (1.12). Then, upon evaluating the effect 
of A, a value of A is selected and used with Eq. (1.13) to evaluate the effect of B. 
This two-step method of evaluating the effect of A and then B may be carried out 
as a two—parameter root locus procedure. First, we obtain a locus of roots as A 
varies, and we select a suitable value of A; the results are satisfactory root locations. 
Then we obtain the root locus for B by noting that the poles of Eq. (1.13) are the 
roots evaluated by the root locus of Eq. (1.15). 

In Mac RootLocus, this design approach is used to calculate the roots for the 


plot. 


Il. OPERATING THE MACROOTLOCUS 


A. BASIC MACINTOSH USE 

This document must serve as a user's manual as well as a technical explanation 
of the program and its capabilities. For this reason, the following brief explanation 
of Macintosh use is included. It is by no means a substitute for the Apple 
Macintosh Users' Manual but it will contain enough information for the beginner to 
be able to use Mac Root Locus. 

1. Basic Operation. 

The Macintosh has a finder, a special application you use to organize and 
manage your document and to start other applications. You use the finder every 
time you start your Macintosh, or whenever you move from one application to 
another. 

The Macintosh screen looks like a light gray desktop, rather than a 
textual list of commands and responses. ‘The desktop simulates the working 
environment. It is initially clean, displaying small graphic images, called icons, 
with short titles directly under them for each disk presently being used and a trash 
can in the lower right corner. An icon is an image representation of an application 
document or a control to a usable function. They offer quick recognition as to the 
type of item they describe and are easier to identify than lists or directories of file 
names with extensions. 

The main interface between the user and Macintosh is the mouse. The 
Macintosh responds instantly to every movement you make with the mouse. You 


can start applications and get documents, work on them, and put them away again 


just by moving the mouse and pressing the mouse button. 

There are three techniques for the mouse: pointing, clicking and dragging. 
Moving the mouse moves the cursor or pointer on the screen in the same direction. 
Positioning the pointer on an item is called pointing to it. The mouse is used to 
select various items or icons on the screen. You select any item or icon to let 
Macintosh know this is what you want to work on next. You select icons by using a 
technique called clicking. As you click the icon, it becomes highlighted. This 
highlighting shows that you select it. This mouse can also be used to drag across 
something. This means the button is pressed and held while the mouse is moved. 
This action is called dragging. When dragging the mouse across the screen, a 
rectangle is outlined. When the button is released, everything within the rectangle 
is now highlighted and selected. Dragging also refers to moving items on the screen. 
This is done by pointing to an item, pressing and holding the button and moving 
the mouse. This will also move an outline of the icon selected and when the button 
is released, the icon 1s moved to the new location. 

Whenever you work with Macintosh, you tell it two things: what you 
want to work on and what you want to do. First, you tell the Macintosh what you 
want to work on by selecting it as you have been doing with icons on the desktop. 
Then you tell the Macintosh what you want to do with the selection. You usually 
do this by choosing a command from a menu. 

Along the top of the screen, in the menu bar, are titles of the menus. 
Pressing the mouse button while you are pointing to a menu title causes the title to 
be highlighted and a menu to appear, much like a window blind being pulled down. 
The menu contains commands you can carry out on what you have selected. 


Commands that you cannot use right now appear dimmed in the menu. When you 


release the mouse button, the menu disappears. 

To choose a command from a menu, you use the same dragging technique 
you used to move icons. As you drag through a menu, each usable command is 
highlighted in turn. If you change your mind about choosing a command, move the 
pointer off the menu and release the mouse button. Nothing is chosen unless you 
release the mouse button while one for the commands is highlighted. You'll follow 
this same pattern whenever you work with the Macintosh; ‘select’ some information, 
the 'choose' an action for it. For example in MacRootLocus, if you double click the 
mouse on MacRootLocus icon, the icon will be selected and MacRootLocus 
application will start running. You can see the menu bar along the top of the screen 
with the greeting message in the center of the screen. Now, you can select and 
choose the items to use MacRootLocus for your design. 

2. Manipulating Window 

The window is the area that displays information on the desk top. You 
view application documents and folders. Folders are a way of storing and organizing 
things, much like in a file drawer. Double clicking on a folder will open a window 
that displays its contents. Folders can be located in disks or within other folders. 
Usually 8 to 12 windows are the maximum that can be open at any one time. When 
several windows are displayed at once, they will usually overlap. If the window you 
want to view is partially covered by another window, pointing anywhere in your 
window and clicking will select it and make it the active window. The active 
window is always in front of all the others. It is the place you want the next action 
to happen, such as move or select icons or open other foiders. 

The active window can be identified by its highlighted title bar with 


narrow horizontal lines on either side of the title. The active window also usually 


has a close box (to close the window) at the top—left corner; at the top—right corner 
is a zoom box that expands the window until it nearly covers the screen. On the 
bottom—right corner is the size box you use to change the size of a window. As you 
click or drag these boxes, you can get the function that you want. 

3. Handling Input 

Most information in the form of data or text is entered through the 
keyboard. The keyboard includes character keys, numerical keys, direction keys and 
other special keys. The return key tells the computer that the data just typed in 
should be accepted now. In MacRootLocus, the apple key in combination with 
another key is often a shortcut to choose a command from a menu. The tap key is 
used to move between data input points in the dialog for entering data. ‘The 
direction key is important in MacRootLocus. It will be used for some word 
processing applications in a small box of the plot window. This will be further 
explained in the next section. 

When the Macintosh requires information from you, it will display a 
dialog box like the one shown in Fig 2.1. It will tell you exactly what data is needed 
and shows you where to enter it. Data insertion points are small boxes in the dialog 
box that let you type in numbers or text. There are usually several such insertion 
points in each dialog box. The tab key lets you move from one point to another to 
enter data. Usually there will already be data in an insertion point box. This is the 
default data. You can change it if you want but you do not have to. When you 
enter new values in an insertion point box, they will become the new default values. 
Data insertion boxes can also be selected: by pointing and clicking with the mouse. 
Clicking the mouse between two characters in the box will let you insert characters 


between them. Double clicking a box will highlight the entire number or an entire 
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Figure 2.1 Sample Dialog 


word if text is entered. You can also select all or portions of numbers or text by 
dragging the mouse across the text you want to select. When all or part of a text is 
selected, it will be highlighted. It can be removed by using the delete or backspace 
key, or it can be replaced by typing in whatever you want. After all data in the 
insertion box is correct, you can enter the data by hitting the return key or by 
clicking the OK button. If you click on Cancel, any changes to the data in the 
insertion boxes will not be saved and the operation which called the dialog box will 


be canceled. If you type in data that does not apply to the insertion box, such as 


10 


typing in letters in the AMin Gain insertion box of one-parameter plot data dialog, 
the Macintosh beeps. It informs you that you inserted the wrong input and you 
should correct the input that you just typed. 

4. Printing Out 

In order to print out your work, there are a few ways available in 
MacRootLocus. One way is to select the print command in the File menu. This 
allows the user to get a hard copy of any plot displayed by MacRootLocus. 

Another way is to do a screen dump which will print the contents of the 
active window immediately. This is done by holding the command key and then 
typing the number '4'. This is a fast way to get print out of a plot in 
MacRootLocus. You can also create a MacPaint document by pressing the 


command and shift key and typing the number '3’. 


You can take up to 10 of these 
'snapshots' and can then alter them with MacPaint for transferring to a word 


processor for lab writeup. 


B. MENUS AND DIALOGS 

There are five menus for MacRootLocus: Apple, File, Edit, Plot and Help. 
Actually, the File, Plot and Help menus are used to get the plot. These have to be 
used in order which will be explained later. 

The Apple and Edit menu are not directly used by MacRootLocus, but they are 
used in order to follow the standard Macintosh programming philosophy for 
user—friendliness. This allows the experienced user to easily adapt to a new 
program since many of the operations are already familiar. Also, this is to allow for 
easy interaction with various other programs and desk accessories [Ref. 1]. 


As mentioned earlier, MacRootLocus presents commands in menus you pull 
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down from the menu bar. As soon as you choose the command you want, the dialog 
box appears for each command. This allows you to insert input data. Now more 
details for menus and dialogs in MacRootLocus will be explained. 
1. Apple menu 

The Apple menu is identified by a small apple in the top left corner. It is 
used primarily for desk accessories. It is also used by most programs as a way of 
offering program information. This is usually the first item of the Apple menu. 

There are several accessories for some applications. For example, the 
Mac's clipboard and scrapbook are used to interact with the word processing 
program. It is highly recommended that the user read the Macintosh Users Guide 
as it explains the use of the Mac's Accessories which will be of use but will not be 
discussed here. 

2. File Menu 

MacRootLocus starts with the file menu first. It offers EQ Parameter, 
Get Coeff, Print Screen, Print Window and Quit commands. These selections also 
have a keyboard shortcut by holding down the command key, which has a clover 
leaf symbol on it, and hitting the first letter of the menu item at the same time. 

a. EQ Parameter 

The EQ Parameter stands for Equation Parameter. As mentioned 
the previous chapter, the MacRootLocus program calculates the roots of the 
characteristic equation for plotting points. It is neccessary to get the degree of 
polynomials and some equation parameters shown in Figure 2.3 since the Laguerre 
algorithm [Ref. 2] is used in order to calculate the roots. 
The degree of the polynomial should be between 1 and 10 in 


MacRootLocus. Then there are the default values for the other parameters. These 


values avoid the convergence error for almost all polynomials. But if convergence 
error messages appear on the screen, you can change these parameter values. Thesc 
parameters must satisfy the following conditions: 

(1) Initial guess > 0 

(2) Maximum Iteration > 100 


(3) Tolerance > 0 


Characteristic Equation Parameter 


Degree of the polynomial | | 


InitGuess 
Maxiter 


Tolerance 


Cancel 





Figure 2.2 Equation Parameter Dialog 


b. Get Coeff 
After getting the degree of the polynomial, you should insert the 
coefficients of the polynomial. The 'Get Coeff’ stands for Get Coefficients. The 


dialog for the coefficients is shown in Figure 2.3. This dialog is varied depending on 
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the degree of the polynomial as shown in Figures 2.3 and 2.4. 

The algebraic expression for the coefficients of the characteristic 
equation of the system may have up to two undetermined parameters (A and B). In 
the case of the one—parameter root locus method, you use only one undetermined 
parameter (A). The routine uses standard algebraic, or infix, notation with 
parenthesis allowed. Operators can include +, —, *, /, and ~ (exponentiation). The 
unary minus sign is allowed. For example, the characteristic equation is 

S3 + (10 + A)S? + (10 * A + 5000 * B)S + 5000 * A =0 


Figure 2.3 shows how you insert these coefficients. 


Characteristic Equation Coefficient Data 


Cancel 


$**2 so] s**0 





Figure 2.3 Characteristic Equation Coefficient Data Dialog 


Box for Third Degree Polynomial 


If you choose the Get Coeff cominand without the degree of the polynomial. 


the message shown in Figure 2.5 appears on the screen. This message tells you that 


degree of the polynomial has not yet been entered. 


Characteristic Equation Coefficient Data 





Figure 2.4 Characteristic Equation Coefficient Data Dialog 


Box for Seven Degree of Polynomial 


There is no Initial Degree. 





Figure 2.5 Message (1) 
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c. Print Screen and Print Window 
These comands allow the user to get a hard copy of any plot 
displayed by MacRootLocus. When you choose the 'Print Screen' command, the 
whole contents of the screen is printed. If you choose the 'Print Window! command, 
the plot that is to be printed should be on the active window of the display. 
Before attempting this method, ensure the printer is properly set up 
for friction feed. 
d.  Quil 
The last item under the File menu is Quit. Selection of this will 
cause you to leave MacRootLocus and return to the desk top. 
3. Edit Menu 
Since these operations are not actually used in MacRootLocus, no further 
explanation of the Edit menu will be presented here. Any additional information 
regarding the Edit menu can be found in the Macintosh Users Manual. 
4. Plot Menu 
There are two commands, One Parameter and ‘Two Parameter. These are 
called to display the dialog boxes of Figures 2.6 and 2.8 to insert the plot data. 
These commands are followed after selecting the EQ Parameter and Get 
Coeff commands. If not, the message shown in Figure 2.7 appears to tell you that 
the degree of polynomial and polynomial coefficients should be entered before you 
choose the plot menu. After reading the message, just click once and this box will 
disappear. 
a. One Parameter 
When you choose the 'One Parameter’ command, the dialog box 


shown in Figure 2.6 for plotting data of the one-parameter root locus method 
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One Parameter Root Locus Plot Data 
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Figure 2.6 One Parameter Root Locus Plot Data Dialog 


appears. The plot default values are shown in Figure 2.6. They can be changed as 
desired. 

First, the user enters the minimum and the maximum gain valucs 
into the 'Min Gain' and the 'Max Gain! insertion box. Next, the user selects one of 
two types of interval, Linear and Logarithmic. When you click the radio button, 
the desired type of interval is chosen. For the 'Linear' interval, the gain step size is 
calculated by subtracting the minimum gain from the maximum gain entered in the 


dialog box, and then dividing by the number of points to plot. Using 'Logarithinic' 


li 


interval can best be described as giving equally spaced intervals on a logarithmic 
scale. 

Most MacRootLocus programs calculate the gain intervals using tlic 
'Linear' interval. This emphasizes gain values that are closer to the max gain. This 
becomes more evident as the maximum gain to the minimum gain ratio increases. 
Using the 'Logarithmic’ interval gives more emphasis to the lower gains so more 
continuous loci can be drawn. As a basic rule of thumb, if the maximum gain to 
minimum gain ratio is greater than 100, selecting 'Logarithmic' interval will give a 


more continuous plot. 


There is no Initial Degree or 


Characteristic Equation Coefficient. 





Figure 2.7. Message Box (2) 


Next, the scale for the axis will be chosen. For the 'Auto Scale’ the 
system calculates the minmum and maximum value of each axis. When the 
‘Manual Scale' is chosen, the user should insert the minimum and maximum valucs 
for each axis. 

The last item, 'Points to Plot' sets the plot resolution. The bigger 
the value you choose, the better resolution plot you get, but the calculation time 


will be longer. 


Two Parameter Root Locus Plot Data 


AMin Gain AMax Gain [250 
BMin Gain a BMax Gain 


@ Linear Point Interval 
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Figure 2.8 Two Parameter Root Locus Plot Data Dialog 


b. Two Parameter 
The 'Two Parameter' command calls the two-parameter plot data 
dialog. It is shown in Figure 2.8. The items shown in Figure 2.8 are similar to 
those shown in Figure 2.6, but several items are different. 
There exists one more undetermined parameter 'B' to be inserted. 
The 'How many loci' item lets you decide how many loci are to be drawn for each 


parameter. The number of loci will be from 1 up to 10 for each parameter. In 
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Figure 2.8, this value is 5. There is no auto—scale for axes. Only the manual scale 
is available. You focus on the interesting area for your design. The last item is tlic 
marking and justification in order to draw the selected 'A' and 'B' values on the 
plot. There are four radio buttons. ‘two buttons are chosen each time for each 
parameter, one for position, the other justification to draw. There are sixteen 
combinations available for this work as shown in Table 2.1. Figure 2.10 shows you 


the ninth case in Table 2.1. 


Table 2.1. The Combination for Marking and Justification 


A Mark neers 
US a 


: 
B Mark 
CS Viel i 
cation 7 L L 


S : Start Point E : Endpoint R: Right HandSide L : Left Hand Side 










QD: Help menu 
MacRootLocus supports an on—line help menu so that the first-time user 
can get desired results without using trial and error. 
Help is the last item in the menu bar. There are the same item names in 
the menu bar. It makes it easy to look for the item for which the user wants 


information. The contents of each item are the subject of section 2.B. 
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6. Information Box 

Finally, MacRootLocus gives you a convenient way to identify your plot. 
It is a small box in which you can type the information you want to memorize. As 
soon as the plot has been completed, the box appears at the bottom of the plot 
window automatically. 

When the up direction key is pressed twice, the cursor comes out on the 
box. Then you can use the keyboard just like a typewriter. The capacity of the box 
is 160 to 200 letters. If you do not want that box, just click once. It will disappear. 


It is shown in Figure 2.10. 
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Sample Plot (1) 


Figure 2.10 
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Il. DETAILED PROCEDURE MODULE DESCRIPTION 


This chapter will basically be programmer's notes covering the significiant 
procedures, functions and libraries to be used in MacRootLocus. MacRootLocus 
follows the Macintosh programming technique to make the program user friendly. 
It is a fairly simple Macintosh application that uses menus, windows, dialog boxcs, 
and graphics. Most procedures and functions called in MacRootLocus were 
developed in Turbo Pascal verson 1.0 for Macintosh. 

MacRootLocus consists of one main program, one resource file and six units. 
The main program integrates the resource file and units, then it shows you the 
menu bar and the greeting message to start the work. The details for these will now 


be explained. 


A. MAIN PROGRAM 
The Main Program is named MacRootLocus.Pas. This program consists of a 
main body and some precedures to handle the systein. 
1. Main Body 
The structure of the main body uses the concept of event driven 


programming. It looks something like this: 


Initialize; set everything up ; 
mepeat keep doing the following 
systemTask; update desk accessories 


if there's an event... 
_,.then handle it 
Nnteeuser 1s done 


if GetNextEvent (everyEvent,theEvent) then 
HandleEvent (theEvent) ; 
until Finished; 
CleanUp; 


tem oe pe pete pete em 


The program is set up once with the user—defined routine 'Initialize.' It 
then enters a loop that continues until some condition (such as the user selecting 
Quit in a menu) causes it to set the boolean flag 'Finished' to true. Within that 
loop, it performs two major tasks. 

First, it calls ‘System Task' (a Toolbox routine), which allows the Mac 
operating system to update any desk accessories that might be in use. Second, it 
calls 'GetNextEvent' (another Toolbox routine) to see if any events have occurred. 
If any have, the highest priority event is returned in the data structure 'the Event.' 
The program then passes the event to ‘Handle Event,’ which is a user—defined 
routine that handles all the different events that might occur. Such events include 
key—presses, selection of menu items, mouse clicks, and windows being opened, 
closed, or resized. When the program is ready to terminate, it calls the 
user—defined routine clean up. 

2. Handle Event Procedure 

When an event occurs, the operating system creates an event record and 
puts it in a queue, ready for you to handle. To see if there is one waiting, you call 
'GetNextEvent,' a boolean function that returns true if there is an event there for 
you. You give it a mask of the events you are interested in; you can use the 
predefined mask 'Every Event! to look at all events. This event is passed to 'Handle 
Event,’ which takes care of it . 'HandleEvent' is just a case statement using the 
'what' field in 'the Event' to determine which of the procedures to call. 

3. DoMouseDown Procedure 
The routine 'DoMouseDown' determines which window the mouse was in when 
the clicking took place and where exactly it happened. Like 'HandleEvent,' 


‘DoMouseDown! is mostly a case statement. 


4. DoUpdate Procedure 

The Macintosh keeps track of a lot of things for you. For one, it tells you 
when some portion of a window needs to be drawn, because of resizing or removing a 
covering window. This is known as an update event, and it requires special 
handling. To handle an update event, the routine 'DoUpdate' saves the current 
'gsrafport' into 'SavePort' and makes 'the Window' the current port so that you can 
write to it. 'BeginUpdate' limits all output to the section of 'theWindow' that 
needs updating. You then do whatever redrawing is needed. When you are done, 
'EndUpdate' lifts those limits, and 'SetPort(SavePort)' restores the old 'grafport’. 

5. | Dokeypress Procedure 

This routine handles the 'Key down' and the 'Autokey' events. It is a 
check to see if a command—key combination was pressed; if so, it checks if the kev is 
a menu command and takes appropriate action. 

6. Handle Menu 

The procedure 'Handle Menu’ decodes the mouse position and figures out 
which menu and which item in that menu were selected. It uses a case statement to 
select the action for the appropriate menu; the menu value is the ID assigned when 
the menu is created. The commands in a menu are numbered from the top down, 
with the first command having a value of one. The action itself is usually a second 
case statement, based on the menu item. 

When an item in a menu is selected, the name of that menu (in the menu 
bar at the top of the screen) is highlighted, that is, inverted to white—on—black. 
When you are done processing the menu command, the menu bar is restored to 
normal by calling 'HiliteMenu(0),' which is at the bottom of ‘HandleMenu.' 


Another procedure is called before you can leave 'HandleMenu,' 'UpdateMenu,' a 


local procedure that tests to see if certain items are to be enabled or disabled. To 
enable and disable menu items, the standard Macintosh procedures 'Enableltem' 
and 'Disableltem' are called. 

7. Initialization 

The initialization procedure for the MacRootLocus program includes the 
following structure: Call 'Init' routines, set up menus, set up windows, do other 
graphic initialization and do program—specific initialization. 

There is an 'Init' routine for most of the major managers. The first, and 
most important, is 'InitGraf' (thePort). That sets up 'QuickDraw!' (which is used 
by just about everything else) and sets up a 'grafport' for the screen. Other 'Init' 
routines are: InitFonts, InitWindow, InitMenus, TEInit, and InitDialogs (NIL). 

Setting up menus involves four steps. First, it defines the menus 
themselves. If a resource file is used, just do a call to 'GetMenu' for each menu 
handle, or even a single call to 'GetNewMBar.' Otherwise, it has to build each 
menu using an initial call to 'NewMenu', followed by a call or calls to Append 
Menu. Second, if it is handling desk accessories, call 'AddResMenu'. Third, add all 
the menus to the menu bar by marking successive calls to 'InsertMenu.' 'Finally, 
call 'DrawMenuBar' to display the menu titles and make them active. 

As with menus, the window initialization takes several steps. If 
MacRootLocus needs a window at start up, create it using either'GetNewwindow' 
(reading in from resources) or Newwindow (building it in place). Having created 
the window, make it the current 'grafport' by calling ‘SetPort', then make it the 
active window by calling ‘Select Window. ' 


The program—specific initialization should probably come here. All of 


the default values for the dialogs are defined and the array vectors are initialized 


here. 


B. GLOBALVAR UNIT 
This unit declares all of the whole variables to be used in MacRootLocus 
except a few of local variables. This unit is called by the main program and all 


units. It defines constant, data type, and variables. 


C. MAKEROOT UNIT 

The unit MakeRoot is very important in MacRootLocus. This unit provides 
several procedures and functions to find the roots and the array vectors for plot. 
since the characteristic equation is derived, the program must be able to parse the 
user's characteristic polynomial coefficient equations in order to understand the 
relations and be able to iteratively substitute in values for undetermined 
parameters: A (for one-parameter root locus method) or A and B (for 
two—parameter root locus method). 

This unit has two main procedures, the Get Root 1 procedure for the 
one—parameter root locus method and the Get Root 2 procedure for two—parameter 
root locus method. The simplified algorithm for the two—parameter root locus 
method is outlined in Figure 3.1. 

There are several procedures to perform this algorithm. These will be 
explained in the following section except for the Rootfinder unit. The InfixtoPolish 
and the ComputePolish are especially interesting. 

1. InfixtoPolish and ComputePolish Procedure 


The coefficients of the polynomial equation are written in algebraic, or 


Le) 
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Given a system's characteristic polynomial : 


CE =a Sep a oO) eae eee 
n n-1 1 0 


where an, 4n-1, etc. are algebraic expressions in A and B. Also A 
is to be stepped from the minimum A to the maximum A and B varied 
from the minimum B to the maximum B and the reverse is processed. 


Then 


SET A = Min A. SET B = Min Bi 
SET A DeltaStep = abs((A Max — A Min) / Step) 
SET B DeltaStep = abs((B Max — B Min) / Step) 
FOR. ket Oued { case of A and B parameter 
FOR. |=" ie torstepado {Step = quantity of losi 
IF Step A then A = A Min + A DeltaStep * (j — 1) 
ELSE B = B Min + B DeltaStep * (j - 1) 
WHILE point no <= (points — 1) do 
IF Step A then B = NextGain2 
ELSE A = NextGain2 
FOR term = Initial Degree downto 0 do 
CONVERT a; from Infix to Polish 
SUBSTITUTE values for A and B 
COMPUTE a; 
END {FOR} 
CALL RootFinder 
CALL Results {make plot array and numerical data} 
END{WHILE} 


A = False {Next case} 
END{FOR} 





Figure 3.1 Two Parameter Root Locus Algorithm 


‘infix' notation. The available operators include (+) addition, (—) subtraction, (*) 
multiplication, (/) division, and (*) exponentiation. These operators follow a 
hierarchical precedence with exponentiation operations being done first, followed by 


multiplication and division, and finally addition and subtraction. Operations like 


multiplication and division which have the same precedence are performed from left 
to right when conflicts arise. To change the order of precedence, parentheses may 
be used around any set of operations. ‘These parenthetical expressions have the 
highest priority and, when nested, the innermost operations within parentheses are 
done first. This scheme follows closely the protocol used in most calculators and 
high level programming languages. 

Infix notation, while convenient for the program user, does not lend itself 
well to computer manipulation. A better way to represent equations for the 
computer is the so called 'reverse Polish notation.' In reverse Polish notation, the 
operands of an equation are entered first, followed by the operator. For example, 
the infix expression 

3* 445 
would be represented as 

oe one 
in reverse Polish notation. The numbers 3 and 4 are entered and multiplied, then 5 
is entered and added to the previous result. Using the concept of a ‘stack’ the 
reverse Polish expression is easy to evaluate. 

Recall that a stack is a last—in—first-out queue whose operation is 
analogous to a stack of trays. ‘To operate the stack, the program calls a 'push' 
procedure to place an item on the stack, and a 'pop' procedure to remove the top 
item. Now, using the example given above, an arithmetic evaluation procedure can 
be illustrated. Figure 3.2 demonstrates such an implementation. 

The basic equation evaluation algorithm can be outlined in three steps: 

(1) Scan the reverse Polish equation term—by—term. 


(2) If the term is a constant then push it onto the stack. 





3 Pomel elle a ee a] > 
<-stack-> {pop4} {pop5} 
{pop3} {pop 12} 
(multiply) L__| {ada} 
{push 3} {push 4} {push 12} {push 17} 


{push 5} 


Figure 3.2. Example of the Stack Operation 


(3) If the term is an operator then pop the first two items off the 
stack, apply the operator, and push the result back onto the 
top of the stack. 

When the algorithm is completed, the answer to the expression will be on the top of 
the stack. 

To get the infix equation into reverse Polish form is a bit more difficult 
than simply evaluating the Polish expression. Of special consideration when 
building the Polish form of the equation are the operator priorities and the use of 
parenthesis to change those priorities. A set of rules can be written which outline 
the conversion procedure [Ref. 2]. These rules are discussed below with an 
illustrative eaxmple. The infix expression 

8+ (7-6/3) *2 


will be parsed using the following operator priority table. 


Operator Priority 
° 4 
one | 2 
a Re Z 
operand 1 
(,) space 0 
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Using a result string called RPN, and an operator stack the rules for infix to reverse 
Polish conversion follow: 

(1) If an operand is encountered, move it to RPN. 

(2) If an operator is encountered, move all higher priority 
operators on the stack to RPN and push the new operator 
onto the stack. 

(3) If a left parenthesis is encountered push it onto the stack. 

(4) If a right parenthesis is encountered, pop all operators off the 
stack and append them to RPN until a left parenthesis is 
encountered. Discard both parenthesis. 

(5) When finished with the infix expression, pop all remaining 
operators from the stack and append them onto RPN. 


Applying these rules to the example problem above is shown in Figure 3.3. 
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Figure 3.3 Conversion from Infix to Reverse Polish 


53 


2. NextGainl Function 
This function makes two kind of intervals, linear interval and 
logarithmic interval. This function is called by the 'GetRoot1' procedure. 
3. NextGain2 Function 
This is like the 'NextGain1l' function except it has two parameters. It 
supports two kind of intervals for each parameter. This function is called by the 
'GetRoot2' procedure. 
4. Results Procedure 
This procedure outputs the calculated roots to the device 'Out File’. 
Then it makes it possible for the user to access the numerical data from the hard 
disk. Also, it supports the plot array named 'GraphArray' for the plotting data. 
5. PlotRootLocus] and PlotRootLocus2 Procedure 
The structure of these procedures is exactly the same except for the 
auto—scale function of the coordinate to be supported by the 'PlotRootLocus1.' 
These procedures draw the roots and some information in the desired coordinate. 
First, the ‘Select Wind' is called with the integer number and a boolean expression in 
order to select a window as visible. The 'Define Header' is called to draw the title. 
Next the 'OpenPic' opens a picture for a specific window and only shows 
the drawing if the boolean is set True. In order to defind the coordinate it calls 
'FindWorld' for the auto—scale or 'FindWorld1l' for the manual—scale. Then 
'DrawAxis' draws an axis with Footers and optional arrows on the axis. Finally, 


'DrawPolygon' draws a polygon defined in the plot array with the predefined shape. 


D. ROOTFINDER UNIT 


The unit which solves for the roots of a polynomial is called 'Root Finder. ' 
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This unit uses Laguerre's method with linear deflation. Since this method is a 
commercially available package of subroutines in the 'Turbo Pascal ToolBox 
(Numerical Methods)’, a brief explanation is offered here. 
1. Laguerre's Method 

Laguerre's method attempts to approximate all the real and complex 
roots of a real or complex polynomial. Laguerre's method is very reliable and quick, 
even when converging to a multiple root. 

To motivate (although not rigorously derive) the Laguerre formulas we 
can note the following relations between the polynomial and its roots and 


derivatives. 


Patt) = aaa —Soe. . (x = oen) (Sa) 


In|/Pa(x)| = Infx — xy] + Infx — xe] +... 


+ In|x — xq| (| 
Gin Ess (o)| 1 2 
x ~ xX — X ae = Xo aes 
l eee 
a ae Se G (335) 


_ _GEloLSdes 5 = se eae) es 


dx? (x - x)? (x -— x9)? 
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Starting from these relations, the Laguerre formulas make what Acton nicely calls 
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"a rather drastic set of assumptions", the root x1 that we seek is assumed to be 
located some distance a from our current guess x, while all other roots are assumed 
to be located at a distance b. 


a=x—xX%; bD=x — x= eee (325) 


Then we can express Eqs (3.3), and (3.4) as 


7+ —=G (3.6) 
eee i (3.7) 
BF b? 


which yields as the solution for a 


gee (3.8) 
Goes tp \(nH — G?) 


where the sign should be taken to yield the largest magnitude for the denominator. 
Since the factor inside the square root can be negative, a can be complex. (A more 
rigorous justification of Eq. (3.8) is in [Ref. 3].) 

The method operates iteratively; for a trial value x, a is calculated by 
Eq. (3.8). Then x — a becomes the next trial value. This continues until a is 
sufficiently small. In the next section, a major procedure which handles this method 


will be explained. 
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2.  InitAndTest Procedure 
This procedure sets the initial value of the input and output variables. 
This procedure also tests the tolerance (Tol), maximum number of iterations 
(MaxIter), and code. Finally, it examines the coefficients of Poly. If the constant 
term is zero, then zero is one of the roots and the polynomial is deflated accordingly. 
Also if the leading coefficient is zero, the degree is reduced until the leading 
coefficient is non—zreo. 
3. FindOneRoot Procedure 
This procedure approximates a single root of the polynomial Poly. The 
root must be approximated within Maxlter iterations to a tolerance of Tol. The 
root, a value of the polynomial at the root, and the number of iterations(Iter) are 
returned. If no root is found, the appropriate error code (Error) is returned. 
4. EvaluatePoly Procedure 
This procedure applies the technique of synthetic division to determine 
the value (yValue), first derivative (yPrime) and second derivative (yDoublePrime) 
of the polynomial, Poly, at X. The Oth element of the first synthetic division is the 
value of Poly at X, the Ist element of the second synthetic division is the first 
derivative of Poly at X, and twice the 2nd element of the third synthetic division is 
the second derivative of Poly at X. 
5. ConstructDifference Procedure 
This procedure computes the difference between approximations; given 
information about the function and its first two derivatives. 
6. ‘TestForRoot Function 
These are the stopping criteria. Four different ones are provided. If you 


wish to change the active criteria, simply comment off the current criteria 
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(including the appropriate OR) and remove the comment brackets from the criteria 
(including the appropriate OR) you wish to be active. 
7. ReducePoly Procedure 
This procedure deflates the polynomial Poly by factoring out the Root. 


Degree is reduced by one. 


E. MESSAGE UNIT 
This unit provides the several messages which inform the user with some 
warnings and help informations. There are 11 procedures,which have the same 
structure, to provide the content of a message and one procedure, 'SetupWindow,' 
to define the window. Here, the 'SetupWindow' procedure and one sample 
procedure to make message will be explained. 
1. SetupWindow Procedure 
This procedure defines all of the windows in the program. It calls the 
'DefineWindow' procedure from 'TurboGraph' Unit. There are five other standard 
window types: documentPrc, DocProc, dBoxProc, PlainDBox, and noGrowDocProc 
[Ref. 4: P. 463]. We can insert window type, window ID, and size of window in this 
procedure. 
2. MakelInfoScreen Procedure 
This is one of nine procedures for messages. It is a sample to explain 
this kinds of procedure. 'MakelInfoScreen' procedure brings up a window vith a 
description of what this program is. First, a 'GrafPort' is called. It is simply a 
Pascal record type with fields that control QuickDraw's behavior.{Ref. 4: p. 405 — 
422| 'GrafPort' allows an application with windows to invoke drawing operations 


that are appropriate for each window. The 'SelectWind' chooses the window to be 
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defined in 'SetupWindow' using just window ID for this procedure. The selected 
window moves to the center of the screen using the 'Movewindow' and the 
‘Set Visibility’ sets the visibility of a window. The 'TextFont,' 'TextSize,' and 
'TextStyle' are used to define the letter to be drawn in the window. Finally the 
'MoveTo' assigns the location to be drawn. 'DrawString' draws the messages in the 


window. 


F. MYDIALOG UNIT 

This unit supports 4 dialog boxes for input data. Each dialog procedure has 
the same skeleton to make a program; the only major procedure to construct the 
program will be explained in this section. 

A dialog box communicates with the background text, controls, icons, and 
'editText' items in which the user can enter and edit text; the user talks back with 
the mouse and keyboard. 

The Dialog Manager [Ref. 5: p. 53] handles this communication. The Dialog 
manager leans heavily on the resource mechanism [Ref. 5: pp. 137 — 154] to provide 
it with data structures. Basically, each dialog is represented on disk as a resource of 
'resType' DLOG—a template describing the dialog's size, window type, and 
title—and a resource of type DITL (Dialog Item List), which lists the content of the 
dialog (controls background text, and so on). 

A program that intends to use dialogs sees to it that appropriate resources of 
type DLOG are available at runtime. Then, the program loads them into memory 
and draws them on the screen with 'GetNewDialog.' 

After 'GetNewDialog' has read a dialog template and its item list into memory 


and drawn it on the screen, a program enters a loop in which it repeatedly 


calls the lynchpin of dialog processing, 'ModalDialog,' and acts on the integer value 
it returns. 'ModalDialog' allows the user to customize its operation with a routine 
specified by the user. The user informs 'ModalDialog' of this by passing a specific 
parameter that points to the routine, which then gets control every time an event is 
generated when the dialog is on the screen. It is up to the filter routine to decide 
what to do about it. 

To set the value of a control, you need to call 'SetCtlValue’. Since the 
'SetCtl Value’ expects a handle to the control record, you will have to get a handle 
to the control that needs setting. This is a task for the Dialog Manager's 
'GetDItem' routine. 'GetDItem' takes in two bits of information and returns three. 
Given the indicated dialog and item number, it returns the information about that 
particular item: its type (such as 'radioButton' or 'staticText', encoded as an 
integer), a handle to its underlying data structure (for controls, this is a 
'ControlHandle'), and finally, its bounding box. The key working with 'editText' 
items is to use 'GetIText' procedures. An 'editText' item begins life with the 
starting content assigned by its definition in an RMaker file [Ref. 5: p. 8]. After the 
user has played with it (as determined by a suitable 'itemHit' value returned by 


'ModalDialog'), 'GetIText' is used to see what the value is now. 


G. TURBOGRAPH UNIT 

This unit supports many procedures for graphics under the Turbo Pascal 
environment. It has commercially available packages in the Turbo Pascal Toolbox 
(Numerical Method). 


Many procedures and functions in this unit are called by several units. Since 
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these routines are already explained, they will not be explained again. The 


comments in source code of this unit give you enough information. 


39 


TV. EXAMPLE OF DESIGN 


A. OVERVIEW 
The root locus method is very valuable in the analysis of dynamic systems, 
and is also used for design. By inspection of the curve, we determine: 
(1) | Whether any loci cross the jw axis into the right half s—plane ( If such a 
crossing exists, it defines a stability limit for the system.), 
(2) The frequency (value of w) at such a stability limit, 
(3) Whether the loci go through any area on the s—plane where we want 
dominant complex roots for our system (i.e., can we get what we want). 
We can also determine the value of the parameter at the stability limit, the 
value of the roots for any specified value of the parameter, and the value of the 
parameter required to place roots at any selected point on the loci. These latter 
items, however, are not done by inspection of the curves, but by inspection of the 
computer printout or by separate calculations. On the MacRootLocus, both plot 
and calculation data are supported. Since two parameters of the system are 
adjustable, it is possible to calculate and plot a family of root loci for the pair of 
parameters The technique for synthesizing a system utilizing the rootlocus method 


is demonstrated in this chapter. 


B. GRAPHICAL SOLUTION 
In this section, four different kinds of control problems will be demonstrated. 


Such problems include simple cascade compensators for systems, subject only to 
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step inputs and/or load disturbances, and feedback compensators with no more than 
two adjustable variables. 
1. Example 1 (cascade lead compensation) 
The block diagram of Figure 4.1 shows the general case of cascade 


compensation. 





Figure 4.1 Simple Casade Compensation Block Diagram 


Consider a plant with 


GA(S)i= Ss) aa (4.1) 
S7(S + 1)(S + 100) 


The uncompensated system is unstable. It is desired to stabilize the system with a 
low—pass filter to reduce bandwidth. Specifications are now in the form of a desired 
location for the roots of the characteristic equation. A cascade lead compensator is 


to be used, and we choose the simplest possible i.e., 


4] 


G.(S) = Bees 4) (4.2) 


Also we define A and B as 
A =—., B = P. (4.3) 


The characteristic equation for this system is 


S° + (101 +B)S* + (100 + 101B)S* + (100B + 400A)S? 
+ 400A(1 + B/A)S + 400B = 0. (4.4) 


The two parameter root locus family for this system is given on Fig 4.2. Since the 
desired roots are 

SO = Nege ee | asi 
an appropriate choice from this plot is A = 5.0, B = 4.8, from which Z = 0.96, and 


P = 5.0. 
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Figure 4.2 Plot of Example 1 


43 


2. Example 2 (cascade lag compensation) 


If we consider a type—one third—order plant for which 


150 
6.(S) = EE (4.5) 
S(S + 1)(S + 10) 
and 
a ( 
€.($) (4.6) 
Z (S + P) 
Then, letting 
hea BS (4.7) 
Z 


the characteristic equation for the cascade lag compensated system is 


St + (11 + B)S? + (10 + 11B)S? + (10B + 150A)S 
+ 150B = 0. (4.8) 


The two parameter root Locus for the lag relocation zone is shown on 
Figure 4.3. Choosing A = 0.1 and B = 0.1 gives P = 0.01, Z = 0.1 and provides 
complex roots at 
a 2 lee 
with real roots at S = —1.014, —10.15. 
Stabilization can also be achieved with a lead filter, but the resulting 


system will have a wider bandwidth and faster response. 
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Figure 4.3. Plot of Example 2 
45 


ay Example 3 (velocity feedback) 
The basic problem is given by the block diagram of Figure 4.4. G(s) 


may be any order. 





Figure 4.4 Velocity Feedback Compensation Block Diagram 


For a second—order system, let 
G(s) = 100 / S(S+2). (4.9) 
We want the root loci for 


EE (4.10) 
S(S + 2) + 100 


so the characteristic equation is 
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The resulting loci are shown on Fig 4.5. Inspection of Fig 4.5 shows that increasing 
the feedback gain increases damping. Any desired value of ¢ is available for the 
complex roots. With high feedback gain, overdamping (all real roots) is available, 
and no positive value of gain can make the system unstable. To design the system, 
pick a root location on the locus, use the magnitude rule to find the gain as in 


Figure 4.5, or, since a computer program was used to generate the locus, the 


tabulated data 


S* + (2 + 100A)S + 100 = 0. 


are supplied to get the desired information. 


data are shown below. 


= 0797 Y 


-98533003622524e+0 
- 98533003622524e+0 


= ORO G 35 


.31655371460862e+0 
.31655371460862e+0 


= O209351 


- 67530563384649e+0 
.67530563384649e+0 


= OaLOlZs 


.06387368606130e+0 
.06387368606130e+0 


= C0969 


-48473591175408e+0 
-48473591175408e+0 


= OF iesesae 


-94057630317454e+0 
-94057630317454et0 


= OB Z369 


-43430192112068et0 
.43430192112068e+0 


+ + 
co 00 CO oO 


+ + 


+ + 
! 
co © 
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.66870719484229e+0 
.66870719484229e+0 


-46960781852863e+0 
-46960781852863e+0 


.23352330186964et0 
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599169390252 /5Z2er0 
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Figure 4.5 Plot of Example 3 
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4. Example 4 (velocity and acceleration feedback) 
The block diagram of Figure 4.6 feedback shows the general case of 


velocity and acceleration feedback. 





Figure 4.6 Velocity and Acceleration Feedback Compensation Block Diagram 


If we consider a third order plant for which 


(4.12) 


then the root locus form is 


2 
S(Aa + Boo} ee: (4.13) 


(S+ 1)? +8 


and the characteristic equation 1s 
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S34 (3 + 8A)S? + (3 + 8B)S+9=0. (4.14) 
The two—parameter root locus family is given on Figure 4.7. We may select a 


desired location for the complex roots such as s = —0.73+ j 1.07, for which A; = 


0.48, By =0.8. The real root is at s = — 5.38 so the complex roots are dominant. 
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Plot of Example 4 
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Figure 4.7 


V. CONCLUSION AND RECOMMENDATION 


MacRootLocus is a useful and powerful tool for the designer and a simple and 
easy to use package for the user. The program supports the two parameter root 
locus method intensively as well as the one parameter root locus method. Also, it 
offers high resolution plots, the tabulated data and a small text editor for recording 
some information of the design. 

The program does, however, have some areas where improvements are 
possible. These are discussed below. 

(1) According to the engineer's needs, the root finder can be replaced 

easily with a more reliable and faster algorithm, and 

i) the capability to move the information window can be added. 

Since this program was developed under the standard interface philosophy the 
Macintosh was designed for, it is a simple matter to change and to append some 
subroutines. 

MacRootLocus is the user friendly CAD program available that is both simple 
enough for the beginning student to easily use as well as powerful and flexible 


enough to benefit the experienced system designer. 


APPENDIX 
SOURCE CODE 


This appendix contains the source code of all modules in MacRootLocus 
except Standard Apple Macintosh libraries. Some subroutines in Numerical Method 
(Turbo Pascal Tool Box) are also included since they are slightly modified. A disk 
is available from Dr. Thaler that contains the MacRootLocus source code and the 


MacRootLocus resource file. 
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MacRootLocus 
version1 .O 


04 DEC 1989 


author 
KO, SUNG HOON 


Il 
a ee a ne a ee ne a ee ne en ee a 


toy 


\S = 
= SYSTEM : Apple Macintosh SE =} 
{= LANGUAGE _ : Turbo Pascal for Macintosh version1.0 =} 
{= LIBRARIES : Numerical Method (Turbo Pascal Tool Box) =} 
— RESOURCE — : MacRootLocus.rsrc =} 
(- =} 


{= This main program handles the system and integrate theresource = 
{= and 6 units. = 


{$B+} { Set the bundle bit } 

{$R MacRootLocus.Rsrc} { Identify resource file for menu and icon info } 
{ST APPLFFTD} { Set the application type and creator } 
{$S+} { Generate segmented code} 

{$I-} { Turn off I/O error checking } 

{$U SpecVar} 

{SU RootsFinder} 

{$U MakeRoot} 

{$U TurboGraph} 

{$U Message} 

{$U MyDialog} 


uses 

MemTypes, QuickDraw, OSIntf, Toollntf,Packintf, 
PasPrinter,SANE,MacPrint, RootsFinder, SpecVar, 

{$S Second Segment} 

TurboGraph, 

Message, 

{$S Third Segment} 

MakeRoot, 

My Dialog; 
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function UpCase(Ch : char) : char; 


inline 

$301F, { UpCase MOVE.W (SP)+,D0 ; GetCh } 
$0C40, 

$0061, { CMP.W #'a’,D0 ; skip if not lower case } 
$6D0A, { BLT.S @1 } 

$0C40, 

$007A, { CMP.W #72’,D0 } 

$6E04, { BGT.S @1 } 

$0440, 

$0020, { SUB.W #$20,D0 } 

$3E80; { @1 MOVE.W DO,(SP) } 


procedure HideAllWindows; 


WindNum : integer; 
begin 
for WindNum := RegendBox downto InfoScrnWind do 
SetVisibility(WindNum, FALSE); 
end; { HideAllWindows } 


procedure ShowAllWindows; 


=sss5555552522552525222==5} 
{= Make the three main windows visible =} 
Meee | = 2 = = = SE ee - -} 
var 

WindNum : integer; 
begin 


for WindNum := RootLocusWind to RegendBox do 
SelectWind(WindNum, TRUE); 
end; { ShowAllWindows } 


procedure MakeWatchCursor; 

var 
CursorHandle : CursHandle; 

begin 
CursorHandle := GetCursor(WatchCursor); 
SetCursor(CursorHandle”’); 

end; { MakeWatchCursor } 


procedure MakeArrowCursor; 
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begin 
InitCursor; 
end; { MakeArrowCursor } 


procedure PrintScreen; 


HideCursor; 
HardCopy(FALSE); 
ShowCursor; 

end; { PrintScreen } 


procedure DoDeskAcc(|Item : Integer); 


var 
SavePort : GrafPtr; 
RefNum — : integer; 
DName _ : string; 


begin 
GeiPort(SavePort); { save port before starting it } 
Getltem(MenuList[AM], ltem, DName); { get name of desk accessory  } 
RefNum := OpenDeskAcc(DName); { and start that sucker up! — }} 
SetPort(SavePort); { restore grafport and continue } 


end; { DoDeskAcc } 


procedure SetltemState(Mndx,Indx : Integer; Flag : Boolean); 


if Flag then 
Enableltem (MenuList[Mndx],Indx) 
else 
Disableltem(MenuList{Mndx],Indx); 
end; { SetltemState } 


procedure HandleMenu(Menulnfo : Longint); 


var 
Menu : Integer; { menu number that was selected _} 
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Item : Inieger; { item in menu that was selected  } 
WindNum : integer; 


begin 
if Menulnfo <> 0 then 
begin 
PenNormal; { set the pen back to normal } 
Menu := HiWord(Menuinfo); { find which menu the command is in } 
Item := LoWord(Menulnfo); { get the command number } 
case Menu of { and carry it out } 
AppiMenu : if Item = 1 then 
DoAbout { bring up “About...” window} 
else 
DoDeskAcc(Iitem); { startdesk accessory} 
FileMenu : case Item of { File menu } 
1 : GetEQParameter; 
2 : GetCoeft; 
3: PrintScreen; {Print screen} 
4: begin 
HideCursor; 
HardCopy(T RUE); { Print top window } 
ShowCursor; 
end; 
5 : Finished := TRUE; { Quit command } 
end; 
EditMenu : case Item of { Edit Menu } 
1..5 : if not SystemEdit(Item-1) then 
{ Do nothing } 
end; 
PlotMenu : case Item of { bring up the dialog } 
1: PlotOneParameter; { for one parameter pilot data} 
2: PlottwoParameter; { for two parameter plot data} 
end; 
HelpMenu : 


case Item of { Help menu } 
1 : InfoGetEQParameter; 
2 : infoGetCoeft; 
3 : infoPlotOneParameter; 
4 : InfoPlotT'woParameter; 
5 : infoPrint 
end 


end;{ case } 
HiliteMenu(0); { reset menu bar } 


Of 


end; 
end; { HandleMenu } 


procedure SelectOurWindow(WPtr : WindowPth); 


PF v4 


== = ea 


begin 
if (WPtr = Window[RootLocusWind].P) then 
SelectWind(RootLocusWind, TRUE); 

end; { SelectOurWindow } 


procedure HandleClick(WPtr : WindowPtr; MLoc : Point); 


begin 
if WPtr <> FrontWindow then { if it’s not in front... } 
SelectOurWindow(WPtr); 
end; { HandleClick } 


procedure HandleGoAway(WPtr : WindowPtr; MLoc : Point); 


{= Handle a mouse click in the go-away box of awindow_ = 
[ess sss sssee eee eee ease ee eee ees} 
var 
WPeek : WindowPeek; { for looking at windows } 
begin 
if WPtr = FrontWindow then { if it’s the active window } 
begin 


WPeek := WindowPeek(WPtr); { peek at the window } 
if TrackGoAway(WPtr, MLoc) then { and the box is clicked } 


begin 
if WPeek*.WindowKind = userKind then { if it's our window } 
HideAllWindows { time to stop } 
else 
CloseDeskAcc(WPeek*.WindowKing) { close DeskAcc } 
end 
end 
else 
SelectOurWindow(WPtr); 


end; { HandleGoAway } 


procedure HandleZoom(WPtr : WindowPtr; MLoc : Point; WLoc : integer); 
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var 
WPeek : WindowPeek; { for looking at windows } 
begin 
if WPtr = FrontWindow then { if it's the active window } 
begin 
WPeek := WindowPeek(WPtr); { peek at the window } 
if TrackBox(WPtr, MLoc, WLoc) then { and the box is clicked } 
begin 
if WPeek*.WindowKind = userKind then { if it's our window } 
ZoomWindow(WPtr, WLoc, FALSE); 
end 
end 
else 
SelectOurWindow(WPtr); 
end; { HandleZoom } 


procedure HandleGrow(WPtr : WindowPtr; MLoc : Point); 


type 
GrowRec = record 
case integer of 
O : (Result : Long!nt); 
1 : (Height, Width : integer); 
end; 
var 
Growlnfo : GrowRec; 
WindNum : integer; 
begin 
if WPtr = FrontWindow then { if it's the active window } 
with Growlnfo do 
begin 
Result := GrowWindow(WPtr, MLoc, GrowArea); { get amount of growth } 
SizeWindow(WPtr, Width, Height, TRUE); { resize window } 
InvalRect(WPtr’.portRect); { set up for update } 
WindNum := 2 ; 
if Window[WindNum].P <> FrontWindow then 
DrawGrowlcon(Window[{WindNum].P); { draw grow icons } 
end 
else 
SelectOurWindow(WPtr); 
end; { HandleGrow } 


procedure HandleDrag(WPtr : WindowPtr; MLoc : Point); 
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var 
WindNum : integer; 
begin 
if WPtr = FrontWindow then 
begin 
DragWindow(WPtr, MLoc, DragArea); { in the drag bar } 
WindNum := 2 ; 
DrawGrowlcon(Window[WindNum].P); { draw grow icons } 
end 
else 
SelectOurWindow(WPtr); 
end; { HandleDrag } 


procedure DoMouseDown(TheEvent : EventRecord); 


var 
theWindow : 
WindowPtr; 
MLoc 
Point; 
WLoc 
integer; 
begin 
MLoc := TheEvent.Where; { get mouse position } 
WLoc := FindWindow(MLoc, theWindow);:{ get the window and the location } 


case WLoc of 
InMenuBar : HandleMenu(MenuSelect(MLoc)); { in the menu } 


InContent : HandleClick(theWindow, MLoc);  { inside the window } 
InZoomiIn : HandleZoom(theWindow, MLoc, WLoc);{ in the zoom box — } 
InZoomOut : HandleZoom(theWindow, MLoc, WLoc);{ in the zoom box} 
InGoAway  : HandleGoAway(theWindow, MLoc); {in the go away box } 
InGrow : HandleGrow(theWindow, MLoc); { in the grow box  } 


InDrag 
HandleDrag(theWindow, MLoc); { in the drag bar  } 


InSysWindow : SystemClick(TheEvent, theWindow); { in a DA window  } 
end; 
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end; { DoMouseDown } 


procedure DoUpdate(TheEvent : EventRecord); 


var 
SavePort, 
theWindow : WindowPtr; 
begin 
theWindow := WindowPtr(TheEvent.Message);  { find which window } 
if (theWindow = Window[RootLocusWingd].P) or 
(theWindow = Window[RegendBox].P) 


then { only update our windows } 
begin 
MakeWatchCursor; 
GetPort(SavePort); { save current grafport } 
SetPort(theWindow); { set aS current port} 
BeginUpdate(theWindow); { signal start of update } 


{ and here's the update stuff! } 


if theWindow = Window[RootLocusWing].P then 

begin 
ClearWindow(RootLocusWind); 
DrawGrowlcon(theWindow); 
DrawPic(RootLocusWind); 

end; 

{ now, back to our program...} 


EndUpdate(theWindow); { signal end of update } 

SetPort(SavePort); { restore grafport } 

MakeArrowCursor; { restore cursor } 
end 


end; { DoUpdate } 


procedure DoActivate(TheEvent : EventRecord); 


var 
AFlag __: boolean; 
theWindow : WindowPtr; 
WindNum _ : integer; 
begin 
with TheEvent do 
begin 
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theWindow := WindowPtr(Message); { get the window } 


AFlag := Odd(Modifiers); { get activate/deactive } 
if AFlag then 
begin { if it's activated... } 

SetPort(theWindow) ; { make itthe port } 
WindNum := RootLocusWind ; 
DrawGrowlcon(Window[WindNum].P); { draw grow icons  } 

end; 

end 


end; { DoActivate } 


procedure DoKeypress(theEvent : EventRecord); 


{= handles keypress (keyDown, autoKey) event =} 
Ba a eS ees eeccossss 
var 

KeyCh : char; 
begin 


KeyCh := Chr(theEvent.Message and charCodeMask); { decode character } 
if (theEvent.modifiers and cmdKey) <> 0 then 
begin { menu key command } 
HandieMenu(MenuKey(KeyCh)) { get menu and item} 
end 
else 
if TextInputEnabled then 
begin 
TEKey(KeyCh, textH); 
TEUpdate(thePort*.portRect, textH); 
end 
else 
sysBeep(1); { do *something” } 
end; { DoKeypress } 


procedure Initialize; 


{= Initialize everything for the program =} 

[== eens ee eee eee e ee see ees] 

var 
Indx 

integer; 
Result : real; 
FileErr : byte; 

begin { initialize all managers used _} 
InitGraf(@thePort); { create a grafport for the screen } 
InitFonts; { start up the font manager } 
InitWindows; { start up the window manager } 


ae 


InitMenus; { start up the menu manager } 
TEInit; 

InitGraphic; 
SetUpWindows; 


TextFont(SystemFont); { initialize the font } 


{ set up menus } 
MenuList[AM] := GetMenu(ApplMenu);{ read menus in from resource fork } 


MenuList{FM] := GetMenu(FileMenu); 
MenuList{EM] := GetMenu(EditMenu); 
MenuList[PM] := GetMenu(PlotMenu); 


MenuList[HM] := GetMenu(HelpMenu); 


AddResMenu(MenuList{AM],,DRVR’'); { pull in all desk accessories } 
for Indx := 1 to 5 do { place menus in menu bar } 
InsertMenu(MenuList[Indx],0); 
for Indx := 1 to 6 do 
SetltemState(EM,Indx,FALSE); { deactivate items in Edit menu  } 
DrawMenuBar; { draw updated menu bar to screen} 
Finished := False; { set program terminator to false } 


{ set drag region } 

SetRect(DragArea, XMinGlb+1, YMinGlb+38, XMaxGlb-1, YMaxGlb-1); 
{ set grow region } 

SetRect(GrowArea, XMinGlb+1, YMinGlb+38, XMaxGlb-1, YMaxGlb-1); 
InitGuess.Re := 1.0: 

InitGuess.Im := 0.0; 

Tolerance := 1e-6; 


Maxlter := 100; 
InitDegree := 0; 
AMinGain := 0.1; 
AMaxGain := 10000; 
BMinGain := 0.1; 
BMaxGain := 10000; 
XMn ‘= -10; 
XMx ‘= 5; 

YMn ‘= -10; 

YMx ‘= 10; 

Step ‘= 4; 

Points := 50; 


FillChar(InitPoly, SizeOf(InitPoly), 0); 
FillChar(xAnswer, SizeOf(xAnswer), 0); 
FillChar(xInitPoly, SizeOf(xInitPoly), 0); 
FillChar(InfixArray, SizeOf(InfixArray), 0); 
MakeWatchCursor; 

MakelnfoScreen; 
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MakeArrowCursor; 
TextlInputEnabled := False; 
InitDegreeStatus := False; 
GetCoeffStatus := False; 


end; { Initialize } 


procedure HandleEvent(TheEvent : EventRecord); 


case TheEvent.What of 
mouseDown 
DoMouseDown(TheEvent); {mouse button pushed  } 
keyDown ~ : DoKeyPress(TheEvent); { key pressed down } 


autoKey : DoKeyPress(TheEvent); { key held down } 

updateEvt : DoUpdate(TheEvent); { window needs updating } 

{activateEvt : DoActivate(TheEvent);} { window made act/inact } 
end 


end; { HandleEvent } 


procedure CleanUp; 


begin 
end; { CleanUp } 


begin { MacRootLocus } 


Initialize; { set everything up } 
repeat { keep doing the following } 
SystemTask; { update desk accessories } 
if GetNextEvent(everyEvent,theEvent) then { _ if there’s an event... } 
HandleEvent(theEvent); { ...then handle it } 
until Finished; { until user is done } 
CleanUp; 


end. { MacRootLocus } 
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unit GlobalVar (5000); 

{= This unit declares the whole variables tobe used = 
{= in MacRootLocus except a part of local variables. = 
{$U RootsFinder} 

interface 

uses 


MemTypes, QuickDraw, OSIntf, Toollntf,Packintf,PasPrinter, 
SANE,MacPrint, RootsFinder; 


const 
InfoScrnWind = 1; { Window number of info window } 
RootLocusWind = 2; { Window number of root locus window } 
RegendBox =3; { Window number of informatiom box window } 
AboutBoxWind = 4; { Window number of about box window } 
AlertBox = 5; { Window number of alert box window } 
HelpWind = 6; { Window number of help window} 
MenuCnt = 5; { total # of menus } 
AppiMenu = 1000; { resource ID of Apple Menu } 
FileMenu = 1001; { resource ID of File Menu } 
EditMenu = 1002; { resource ID of Edit Menu } 
PlotMenu = 1003; { resource ID of Plot Menu } 
HelpMenu = 1004; { resource ID of Help Menu } 
AM = 1;  { index into MenuList for Apple Menu } 
FM = 2; {ditto for File Menu } 
EM = 3; {ditto for Edit Menu } 
PM = 4; {ditto for Plot Menu } 
HM = 5; {ditto for Help Menu } 
MaxPlotGlb = 1024;{ The maximum number of points in a plot array } 
type 


str255 = string [255] ; 

StringArray = array[1..20] of str255; { Storage for plot data } 

NodePtr = “Node; 

Node = 
record 
Data : Extended; 
Next : NodePtr; 
end; 

PlotArray = array[1..MaxPlotGlb, 1..2] of Extended; 
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var 
Answer, polish : str255; 
InfixArray : StringArray; 
A,B,ADeltaStep, Alncrem ,BDeltaStep, Blncrem: Extended; 


InitGuess : TNComplex; { Initial approximation } 
Tolerance : Extended; { Tolerance of approximation } 
AMaxGain, BMaxGain : Extended; { Range of unknown parameter values } 


AMinGain, BMinGain : Extended; 
Root, Imag, Value, Deriv : TNvector; { Resulting roots and other info } 


lter : TNintVector; { Iterations to find each root } 

Maxiter : integer: { Maximum number of iterations } 
Initbegree, Degree : integer; { Initial and final degree } 

xInitPoly : TNvector; { of polynomial } 

InitPoly, Poly : TNCompVector; { Initial and final coefficients } 


{ of the polynomial } 
xAnswer : TNCompVector; 
yAnswer : TNCompVector; 
NumRoots, Arraylndex, pointno : integer; { Number of roots } 
Error : byte; { Error flag } 
StepA, StepB, linear : Boolean; 
Step,Points : Longint; 


OutFile — : text; { The output file used by an application} 
OutName — : string; { The out file name} 

lOerr : boolean; { Flags I/O errors } 

Finished _ : boolean; { used to terminate the program } 
theEvent : EventRecord; { event passed from operating system } 


MenuList — : array[1..MenuCnt] of MenuHandle; { holds menu info } 
DragArea : Rect ; { Area in which window can be dragged } 
GrowArea : Rect; { Area in which a window’s size can change } 


DataPicture : PicHandle; { Holds a picture of the plotted data } 
theDialog : DialogPtr: 

itemHit  : Integer; 

theType _ : Integer; 

r : Rect; 

done, AutoScale _—_: Boolean; 
n : Integer; 
h,h3,h4,n5,he : Handle; 
S > Str255; 

XMn : Extended; 

XMx : Extended; 

YMn : Extended; 

YMx : Extended; 


GraphArray : “PlotArray; 
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Ds : DecStr ; 
ARJustification, AMarkStatus,BRJustification, BMarkStatus : Boolean; 


Da, Db, Dan, Dax, Don, Dbx, Dss : DecStr: 
txRect : Rect; 

textH : TEHandle; 

TextinputEnabled : Boolean; 
InitDegreeStatus : Boolean; 
GetCoeffStatus : Boolean; 


implementation 


begin 
end.{GlobalV ar} 
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unit MakeRoot(3000); 


{= This unit provides several procedures and functions to 

{= make the roots and the array vectors for plot. 

SSS SS SS SS SS a Se aS eS Se See Se SS SSeS = 
{$l-} { Disabie I/O error trapping } 

{$S+} { Enable segmentation of code } 


{$U SpecVar} 
{$U Message} 
{$U RootsFinder} 
{$U TurboGraph} 


interface 


uses 
MemTypes, QuickDraw, OSIntf, Toollntf, Packlntf,PasPrinter, 
SANE, MacPrint,RootsFinder, SpecVar, 


{$S SecondSegment} 
TurboGraph, 
Message; 


var 
Error : byte; 


procedure GetRoott1; 
procedure GetRoot2; 


implementation 
procedure InfixToPolish( Answer:str255;var RPN:str255); 


i This procedure converts an algebraic expression(infix) 
{= to reverse Polish notation. 


var 
Stack =: Array[1..50] of char; 
Top,I,N,p : Integer; 

Ch : char; 

Chit : String; 
Firstchar,PrevDigit : Boolean; 


Function Priority(Ch:char): Integer; 
Begin 
Case Ch of 


: Priority:= 4; 


: Priority:= 3; 
1 : Priority:= 3; 
a : Priority:= 2; 
‘+! : Priority:= 2; 


‘a’..'b' : Priority:= 1; 
‘A’..'B' : Priority:= 1; 


'0’..'9' : Priority:= 1; 


% : Priority:= 1; 
ic : Priority:= 0; 
a : Priority:= 0; 
- : Priority:= 0; 
end; 

end; 

begin 

RPN := "; 

io —0; 


FirstChar :=True; 
PrevDigit := False; 


for i:= 1 to Length(Answer) do 


begin 
chi := copy(Answer,i,1); 
ch:= chi[1]; 


P:= Priority(ch); 
if Firstchar and (Ch = ‘-') then P:= 1; 


ifP = 1 then 
begin 
if PrevDigit then 
RPN := concat(RPN,ch) 
else 
RPN := concat(RPN,' ',ch); 


firstchar := False; PrevDigit := True; 
end; 


if P>1 then 
begin 


while (Top>0) and (Priority(stack[top]) >= P) do 
begin 


RPN := concat(RPN,’ ',Stack[Top)}); 
Top := Top - 1; 
end; 
Top := Top + 1; 
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Stack[Top] := ch; 
Firstchar := True; PrevDigit := False; 


end; 
ifch =’(’ then 
begin 

Top := Top +1; 


Stack[Top] := ’(’; 
FirstChar := True; PrevDigit :=False; 
end; 
if ch =’)’ then 
begin 
while Stack[Top]<> ‘(’ do 
begin 
RPN := concat(RPN,’ ’,Stack[Top)); 
Top := Top -1; 
end; 
Top := Top -1; 
FirstChar := False; PrevDigit := False; 
end; 
end; 
while Top > 0 do 
begin 
RPN := concat(RPN,’ ',Stack[Top}); 
Top := Top -1; 
end; 
end; 


procedure ComputePolish(polish:str255;a,b:Extended; var EvalArray:Extended); 


{= This procedure uses the string generated in Procedure =} 
{= InfixToPolish to evaluate the numeric expression =} 
= S66 S655 946555 55s5565 =5=55555555= 5) 
var 

i,N : integer; 

ch,ch2 : char; 

ch1 : String; 

temp : String[255}; 


Value1,Value2,Value3 : Extended; 
StackPtr : NodePtr; 
StackEmpty : Boolean; 


procedure CreateStack; 


New(StackPtr); 
with StackPtr* do 


begin 

Next := nil; 

Data := 0.0; 
end; 

StackEmpty := True; 
end; 


procedure Pop(var Val : Extended);{Push a number onto numeric stack} 


var 
NPtr :NodePtr; 
begin 
if not StackEmpty then 
begin 
NPtr := StackPtr’.Next; 
StackPtr*.Next := NPtr*.Next; 
Val := NPtr*.Data; 
Dispose(NPtr); 
StackEmpty := (StackPtr’*.Next = nil); 
end; 
end; 


procedure Push(Val : Extended); 


{= pop a number off numeric stack a 
ee es oe ee 2 2 = SS} 
var 
NPtr : NodePtr; 
begin 
StackEmpty := False; 
New(NPtr): 
NPtr’.Data := Val; et 


NPtr’.Next := StackPtr’.Next; 
StackPtr’ Next := NPtr; 


end; 


procedure DeleteStack;{Delete stack} 


var 
Temp : Extended; 


begin 
while not StackEmpty do 
Pop(Temp); 
Dispose(StackPtr); 
end; 


function Expon(y,x:Extended) :Extended; 


{= Computes Y raised to X power = 
Bano SSeS so eS soe ss54] 
begin 
Expon := exp( x * (In(y))); 
end; 
begin 
CreateStack; {Initialize} 
temp :=''; 


for i:= 1 to Length(polish) do{do one char at a time} 
begin 

chi := copy(polish,i,1); 

ch:= chi[1];{Get a char} 


case ch of fand evaluate it} 
'0'..'9° : temp :=concat(temp,ch);{Real constant} 
s : temp ‘= concat(temp,ch); 
: begin 
ch1 := copy(polish,i+1,1); 
ch2:= ch1[{1];{Get a char} 


if (ch2 <> '’) and (i < length(polish)) then{Unary minus} 
temp := concat(temp,ch) 
else 
begin {minus operator} 
POP(Value1);POP(Value2); 
Value3 := Value2 - Value; 


PUSH(Value3); 
end; 
end; 
‘a’ : PUSH(a); 
A’ : PUSH(A); 
'b’ : PUSH(b); 
‘iB >: PUSH(B); 


(2 


i 


end; 
end; 


‘begin 
POP(Value1); POP(Value2); 
Value3 := Value2 + Value; 


PUSH(Value3); 
end; 


‘begin 
POP(Value1); POP(Value2); 
Value3 := Value2 * Value1; 
PUSH(Value3); 


end; 


‘begin 
POP(Value1); POP(Value2); 
Value3 := Value2 / Value; 
PUSH(Value3); 

end; 


‘begin 
POP(Value1); POP(Value2); 
Value3 :=Expon( Value2 , Value1); 
PUSH(Value3); 

end; 


begin 
if temp <> ' ‘then 
begin 
Value1 := Str2Num(temp); 
PUSH(Value1); 
temp :=''; 
end; 
end; 


if temp <> '' then 


begin 
Value3 := Str2Num(temp); 


PUSH(Value3): 


end; 


POP(Value1); 


as 


EvalArray := Value; 
DeleteStack; 
end; 


function Ten2 (power: Extended) : Extended; 
begin 

Ten2 := Exp(power * Ln(10)); 
end; 


function NextGain1 : Extended; 


} 
= This function makes two kind of intervals for one =} 
{= parameter root locus method. } 


var 
GainOut, logmin, logmax : Extended; 


begin 
if linear then 
begin 
Alncrem := abs((AMaxGain - AMinGain)/(points - 1)); 
GainOut := AMinGain + Alncrem * pointno; 


end 
else 
begin 
if (AMinGain < 1E-5) and (AMaxGain > 1E-1) then 
logmin := -5 
else 


logmin := Ln(AMinGain)/Ln(10); 
logmax := In(AMaxGain)/Ln(10); 
Alncrem := (logmax - logmin)/(Points -1); 
GainOut := Ten2(logmin + pointno * Alncrem); 
end; 
NextGain1 := GainOut; 
end; 


{= This function makes two kind of intervals for each = 
le parameter of two parameter root locus method. = 


var 
GainOut, logmin, logmax : Extended; 


begin 


if StepA then 
begin 
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if linea: then 
begin 
Bincrem := abs((BMaxGain - BMinGain)/(points - 1)); 
GainOut := BMinGain + Bincrem * pointno; 


end 
else 
begin 
if (BMinGain < 1E-5) and (BMaxGain > 1E-1) then 
logmin := -5 
else 
logmin := Ln(BMinGain)/Ln(10); 


logmax := In(BMaxGain)/Ln(10); 
Bincrem := (logmax - logmin)/(Points -1); 
GainOut := Ten2(logmin + pointno * Bincrem); 
end; 
NextGain2 := GainOut; 
end 
else 
begin 
if linear then 
begin 
Alncrem := abs((AMaxGain - AMinGain)/(points - 1)); 
GainOut := AMinGain + Alncrem * pointno; 


end 
else 
begin 
if (AMinGain < 1E-5) and (AMaxGain > 1E-1) then 
logmin := -5 
else 


logmin := Ln(AMinGain)/Ln(10); 

logmax := In(AMaxGain)/Ln(10); 

Alncrem := (logmax - logmin)/(Points -1); 
GainOut := Ten2(logmin + pointno * Alncrem); 


end; 
NextGain2 := GainOut: 
end; 
end; 
procedure Results(A,B : Extended; 


xAnswer : TNCompVector; 
Error : byte); 


This procedure outputs the results to the device OutFile = 
{= and make the array for plot data. = 


var 
Term : integer; 


begin 
if Degree > 0 then 
begin 
Writeln(OutFile,'The deflated polynomial:’); 
for Term := Degree downto 0 do 
Writeln(OutFile, ‘Poly[',Term,'):’, 
Poly[Term].Re:23, ’ +’, Poly[Term].lm:23,'7’); 
end; 


if Error <= 1 then 
begin 
Writeln(OutFile); 
Writeln(OutFile,’ A =’, A :10:5,'B =', B :10:5 ); 
for Term := 1 to NumRoots do 
begin 
Writeln(OutFile, xAnswer[Term].Re:23, ’ +’, xAnswer[Term].lm:23, ’I'); 
GraphArray“[Array|lndex,1] := xAnswer[Term].Re; 
GraphArray“[Arraylndex,2] := xAnswer[Term].Im; 
Arrayindex := Arraylndex + 1; 
end; 
end; 
end; { procedure Results } 


procedure PlotRootLocus1; 


begin 
SelectWind(2, TRUE); 
DefineHeader(2, 'One Parameter RootLocus’); 
DrawGrowlcon(Window/[2].P); 
OpenPic(2, TRUE); { Open up a picture for window #1 } 


{ Draw the axis and data points } 
if AutoScale then 
FindWorld(2, GraphArray*, Arraylndex) 
else 
FindWorld1(2,XMn, YMn,XMx, YMx); 
DrawAxis( Da , ", True); 
DrawPolygon(GraphArray’, 1, -Arraylndex, -3, 1, 0, TRUE); 
ClosePicture; { Close the picture for window #2} 
ValidRect(Window[2].P*.portRect); 
end; { PlotPowerExpLog } 


procedure PlotRootLocus2; 


begin 
SelectWind(2, TRUE); 
DefineHeader(2, ‘Two Parameter RootLocus’); 
DrawGrowlcon(Window[2].P); 
OpenPic(2, TRUE); { Open up a picture for window #2 } 


{ Draw the axis and data points } 
FindWorld1(2,XMn,YMn,XMx, YMx); 
DrawAxis(Da, Db, True); 
DrawPolygon(GraphArray*, 1, -Arraylndex, -3, 1, 0, TRUE); 
ClosePicture; { Close the picture for window #2 } 
ValidRect(Window([2].P*.portRect); 
end; { } 


procedure GetRoot1; 


{= This procedure gets the root of the characteristic equation for =} 
— one parameter root locus method. =} 
a le = = = 
var 
term : integer; 
f : DecForm; 
begin 


Num2Str(f,AMinGain,Dan); 
Num2Str(f,AMaxGain,Dax); 
Da := ’A:’+ Dan +’ To ’+Dax; 

Ds :='': 
Arraylndex := 1; 
New (GraphArray); 


pointno := 0; 
while pointno <= ( Points - 1) do 
begin 
A := NextGain1; 
B := 0; 
for Term := InitDegree downto 0 do 
begin 


infix ToPolish(InfixArray[Term],polish); 
ComputePolish(polish,a,b,xInitPoly[Term)}); 
InitPoly[Term].Re := xInitPoly[Term]; 
InitPoly[Term].lm := 0.0; 


end; 
Degree := InitDegree; 
Poly := InitPoly; 
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Laguerre(Degree, Poly, InitGuess, Tolerance, Maxlter, 
NumRoots, xAnswer, yAnswer, Iter, Error); 
Results(A, B, xAnswer, Error); 
pointno := pointno + 1; 
end; 
Arraylndex := Arraylndex - 1; 
PlotRootLocus1; 
end; 


procedure GetRoot2{(var StepA : boolean; A,B : Extendeda)}; 


{= This procedure gets the root of the characteristic equation for =} 
= two parameter root locus method. = 
{2-2 22555555525555555555 e555 s5555555]) 
var 
i, J, K, term : integer; 
f : DecForm; 
begin 


ADeltaStep := abs((AMaxGain - AMinGain)/(Step)); 
BDeltaStep := abs((BMaxGain - BMinGain)/(Step)); 
Num2Str(f,AMinGain,Dan); 
Num2Str(f,AMaxGain,Dax); 

Da := ’A:'+ Dan +’ To ’+Dax; 

Num2Str(f, BMinGain,Dbn); 
Num2Str(f,BMaxGain,Dbx); 

Dod := ’B :’+ Dbn +’ To '+Dbx; 

New (GraphArray); 


StepA := True; 
fori := 1 to 2 do 
begin 
for j := 1 to Step do 
begin 


Array|lndex := 1; 
if StepA then A := AMinGain + ADeltaStep*(j-1) 
else B := BMinGain + BDeltaStep*(j-1); 


if StepA then 
begin 
Mum2Sir(f,A,Dss); 
Ds := ’A=’+ Dss; 
end 
else 
begin 
Num2Str(f,B,Dss); 
Ds := 'B=’+ Dss; 
end; 
pointno := 0; 


while pointno <= ( Points - 1) do 
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begin 
if StepA then B := NextGain2 
else A := NextGain2; 
for Term := InitDegree downto 0 do 
begin 
infixToPolish(InfixArray[Term],polish); 
ComputePolish(polish,a,b,xInitPoly[T erm)); 
InitPoly{Term].Re := xInitPoly[Term]; 
InitPoly[Term].lm := 0.0; 


end; 
Degree := InitDegree; 
Poly := InitPoly; 


Laguerre(Degree, Poly, InitGuess, Tolerance, Maxlter, 
NumRoots, xAnswer, yAnswer, Iter, Error); 
Results(A, B, xAnswer, Error); 
pointno := pointno + 1; 
end; 
Arraylndex := Arraylndex - 1; 
PlotRootLocus2; 
end; 
StepA := False; 
end; 
end; 
begin 
end. { Unit MakeRoot } 


unit RootsFinder(2000); 


{= Turbo Pascal Numerical Methods Toolbox =} 
{= Copyright (C) 1987 Borland International } 
is z, 
{= This unit provides procedures for finding the roots =} 
(- of a single equation in one real variable. =} 

} 


{$R+} { Enable range checking } 
interface 


uses 
MemTypes; 


const 
TNNearlyZero = 1E-015; { Close to zero } 
TNArraySize = 10; { Maximum size of vectors } 


type 
TNvector = array(0..TNArraySize] of Extended; 


TNintVector = array[0..TNArraySize] of integer; 


TNcomplex = record 
Re, Im : Extended; 
end; 


TNCompVector = array[0..TNArraySize] of TNcomplex; 


procedure Laguerre(var Degree _ : integer; 
var Poly  :TNCompVector; 
Initauess : TNcomplex; 
Tol ; Extended; 
Maxlter : integer; 
var NumRoots : integer; 
var Roots : TNCompVector; 
var yRoots :  TNCompVector; 
var Iter : TNintVector; 
var Error _ : byte); 


= Input: Degree, Poly, Init@uess, Tol, Maxlter = 
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Qutput: Degree, Poly, NumRoots, Roots, yRoots, Iter, Error 


Purpose: This unit provides a procedure for finding all the 
roots (real and complex) to a polynomial. 
Laguerre's method with deflation is used. 
The user must input the coefficients of the quadratic 
and the tolerance in the answers generated. 


Pre-defined Types: TNcomplex = record 
Re, Im : Extended; 
end; 


TNintVector = array[0..TNArraySize] of integer; 
TNCompVector = array[0..TNArraySize] of TNcomplex; 


Variables: Degree _ : integer; degree of deflated 
polynomial 
Poly |: TNCompVector; coefficients of deflated 
polynomial where Poly[n] is 
the coefficient of Xn 
Initauess : TNcomplex; _ initial guess to a root 
(may be very crude) 
Tol : Extended; _ tolerance in the answer 
Maxlter : integer; number of iterations 
NumRoots : integer; number of roots calculated 
Roots : TNCompVector; the roots calculated 
yRoots :TNCompVector; the value of the function 
at the calculated roots 


Iter - TNintVector; number iteration it took to 
find each root 
Error _: byte; flags an error 


Errors: 0: No error 
2: Degree <= 0 
3: Tol <= 0 
4: Maxlter < 0 


The following inline procedure and function are used to call the user 


defined procedures and functions pointed to by the ProcAddr parameter. 


ol 
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} 
=} 


function UserFunction(X : Extended; ProcAddr : ProcPtr) : Extended; 
inline 

$205F, { MOVE.L (A7)+, AO } 

$4E90; {JSR (AO) } 


procedure UserProcedure(X : TNcomplex; var Y : TNcomplex; ProcAddr : ProcPtr); 
inline 

$205F, { MOVE.L (A7)+, AO } 

$4E90; {JSR (AO) } 


procedure Laguerre{(var Degree : integer; 
var Poly : TNCompVector; 
InitGuess : TNcomplex; 
Tol : Extended; 
Maxliter : integer; 
var NumRoots : integer; 
var Roots : TNCompVector; 
var yRoots :  TNCompVector; 
var Iter : TNintVector; 
var Error _ : byte)}; 


type 
TNquadratic = record 
A, B, C : Extended; 
end; 


var 
Addlter _: integer; 
InitDegree : integer; 
InitPoly : TNCompVector; 
GuessRoot : TNcomplex; 


{----------- Here are a few complex operations. ------------ } 


procedure Conjugate(var C1, C2 : TNcomplex); 
begin 

C2.Re := C1.Re; 

C2.Ilm := -1.0 * C1.Im; 
end; { procedure Conjugate } 


function Modulus(var C1 : TNcomplex) : Extended; 
begin 

Modulus := Sqrt(Sqr(C1.Re) + Sqr(C1.Im)); 
end; { function Modulus } 


procedure Add(var C1, C2, C3 : TNcomplex); 
begin 

C3.Re := C1.Re + C2.Re; 

C3.Im := Ci.lm + C2.1m; 
end; { procedure Add } 


procedure Sub(var C1, C2, C3 : TNcomplex); 
begin 

C3.Re := C1.Re - C2.Re; 

C3.Im := C1.lm - C2.Im; 
end; { procedure Sub } 


procedure Mult(var C1, C2, C3 : TNcomplex); 
begin 
C3.Re := C1.Re * C2.Re - C1.lm * C2.1m; 
C3.lm := C1.lm * C2.Re + C1.Re * C2.1m; 
end; { procedure Mult } 


procedure Divide(var C1, C2, C3 : TNcomplex); 
var 
Dum1, Dum2 : TNcomplex; 
E : Extended; 
begin 
Conjugate(C2, Dum1); 
Mult(C1, Dum1, Dum2); 
E := Saqr(Modulus(C2)); 
C3.Re := Dum2.Re / E; 
C3.lm := Dum2.im / E; 
end; { procedure Divide } 


procedure SquareRoot(var C1, C2 : TNcomplex); 
var 
R, Theta : Extended; 
begin 
R := Sart(Sqr(C1.Re) + Sqr(C1.lm)); 
if ABS(C1.Re) < TNNearlyZero then 
begin 
if C1.lm <0 then 
Theta := Pi/2 
else 
Theta := (-1.0 * Pi) / 2; 
end 
else 
if C1.Re < 0 then 
Theta := ArcTan(C1.lm / C1.Re) + Pi 
else 
Theta := ArcTan(C1.lm / C1.Re); 
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C2.Re := Sart(R) * Cos(Theta / 2); 
C2.Im := Sart(R) * Sin(Theta / 2); 
end; { procedure SquareRoot } 


procedure InitAndTest(var Degree : Integer; 
var Poly : TNCompVector; 
Tol : Extended; 


Maxiter : integer; 
InitGuess : TNcomplex; 
var NumRoots_: integer; 
var Roots  : TNCompVector; 
var yRoots :TNCompVector; 
var Iter : TNintVector; 
var GuessRoot : TNcomplex; 
var InitDegree : integer; 
var InitPoly : TNCompVector; 
var Error : byte); 


{= Input: Degree, Poly, Tol, Maxiter, InitGuess 
= Output: InitDegree, InitPoly, Degree, Poly, NumRoots, 
= Roots, yRoots, Iter, GuessRoot, Error 


{ 
{ This procedure sets the initial value of the above 

{= variables. This procedure also tests the tolerance 

{= (Tol), maximum number of iterations (Maxlter), and 
{= code. Finally, it examines the coefficients of Poly. 

{ If the constant term is zero, then zero is one of the 


{= roots and the polynomial is deflated accordingly. Also 

— if the leading coefficient is zero, then Degree is 

{= reduced until the leading coefficient is non-zero. 
{eeezeeezeeee2eeeeeee2zeee2eee2e2e 222222222252 
var 


Term : integer; 


begin 
Error := 0; 
if Degree <= 0 then 
Error := 2; { degree is less than 2 } 
if Tol <= 0 then 
Error := 3; 
if Maxiter < 0 then 
Error := 4; 


if Error = 0 then 
begin 


$4 
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NumRoots := 0; 
GuessRoot := InitGuess; 
InitDegree := Degree; 
InitPoly := Poly; 
{ Reduce degree until leading coefficient <> zero } 
while (Degree > 0) and (Modulus(Poly[Degree]) < TNNearlyZero) do 
Degree := Pred(Degree); 
{ Deflate polynomial until the constant term <> zero } 
while (Modulus(Poly[0}]) = 0) and (Degree > 0) do 
begin 
{ Zero is a root } 
NumRoots := Succ(NumRoots); 
Roots[NumRoots].Re := 0; 
Roots[NumRoots].lm := 0; 
yRoots[NumRoots].Re := 0; 
yRoots[NumRoots].Im := 0; 
lter[NumRoots] := 0; 
Degree := Pred(Degree); 
for Term := 0 to Degree do 
Poly[Term] := Poly[Term + 1]; 
end; 
end; 
end; { procedure InitAndTest } 


procedure FindOneRoot(Degree _ : integer; 
Poly :  TNCompVector; 
GuessRoot : TNcomplex; 
Tol : Extended; 
Maxlter : integer; 
var Root  : TNcomplex; 
var yValue : TNcomplex; 
var Iter : integer; 
var Error — : byte); 


a = 
{= Input: Degree, Poly, GuessRoot, Tol, Maxlter = 
= Output: Root, yValue, Iter, Error » 


am, 
ll 
II 


{= This procedure approximates a single root of the polynomial 
{= Poly. The root must be approximated within Maxlter 

= iterations to a tolerance of Tol. The root, value of the 

{= polynomial at the root (yValue), and the number of iterations 
=  (lter) are returned. If no root is found, the appropriate error 
{= code (Error) is returned. 


Reet Nee Nee Need ee eee es tee ey te 
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Found : boolean; 
Dif : TNcomplex; 
yPrime, yDoublePrime : TNcomplex; 


procedure EvaluatePoly(Degree : Integer; 
Poly : TNCompVector; 
X : TNcomplex; 
var yValue : TNcomplex; 
var yPrime : TNcomplex; 
var yDoublePrime : TNcomplex); 


{= Input: Degree, Poly, X 

{= Output: yValue, yPrime, yDoublePrime 

{= This procedure applies the technique of synthetic division to 

{= determine value (yValue), first derivative (yPrime) and second 
{= derivative (yDoublePrime) of the polynomial, Poly, at X. 

{= The Oth element of the first synthetic division is the 

{= value of Poly at X, the 1st element of the second synthetic 

{ division is the first derivative of Poly at X, and twice the 

{= 2nd element of the third synthetic division is the second 

{= derivative of Poly at X. 

{ 


var 
Loop : integer; 
Dummy, yOPdummy : TNcomplex; 
Deriv, Deriv2 : TNCompVector; 


begin 
Deriv[Degree] := Poly[Degree]; 
for Loop := Degree - 1 downto 0 do 
begin 
Mult(Deriv[Loop + 1], X, Dummy); 
Add(Dummy, Poly[Loop], Deriv[Loop]); 
end; 
yValue := Deriv[O]; { Value of Poly at X } 


Deriv2([Degree] := Deriv[Degree]; 
for Loop := Degree - 1 downto 1 do 
begin 
Mult(Deriv2[Loop + 1], X, Dummy); 
Add(Dummy, Deriv[Loop], Deriv2[Loop)}); 
end; 
yPrime := Deriv2[1]; { 1st deriv. of Poly at X } 


yDPdummy := Deriv2(Degree]; 
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for Loop := Degree - 1 downto 2 do 
begin 
Mult(yOPdummy, X, Dummy); 
Add(Dummy, Deriv2[Loop], yOPdummy); 
end; 
yDoublePrime.Re := 2 * yDPdummy.Re; { 2nd derivative of Poly at X } 
yDoublePrime.lm := 2 * yOPdummy.|m; 
end; { procedure EvaluatePoly } 


procedure ConstructDifference(Degree : integer; 
yValue : TNcomplex; 
yPrime : TNcomplex; 
yDoublePrime : TNcomplex; 
var Dif : TNcomplex); 


{= Input: Degree, yValue, yPrime, yDoublePrime 
= Output: Dif 


{= This procedure computes the difference between approximations; 
= given information about the function and its first two 
{= derivatives. 


var 
yPrimeSQR, yTimesyDPrime, Sum, SRoot, 
Numert, Numer2, Numer, Denom : TNcomplex; 


begin 

Mult(yPrime, yPrime, yPrimeSQR); 
yPrimeSQR.Re := Sqr(Degree - 1) * yPrimeSQR.Re; 
yPrimeSQR.Im := Sqr(Degree - 1) * yPrimeSQR.Im; 

Mult(yValue, yDoublePrime, yTimesyDPrime); 
yTimesyDPrime.Re := (Degree - 1) * Degree * yTimesyDPrime.Re; 
yTimesyDPrime.Im := (Degree - 1) * Degree * yTimesyDPrime.|m; 
Sub(yPrimeSQR, yTimesyDPrime, Sum); 
SquareRoot(Sum, SRoot); 

Add(yPrime, SRoot, Numer1); 

Sub(yPrime, SRoot, Numer2); 

if Modulus(Numer1) > Modulus(Numer2) then 

Numer := Numer 
else 

Numer := Numer2; 
Denom.Re := Degree * yValue.Re; 
Denom.Im := Degree * yValue.Im; 

if Modulus(Numer) < TNNearlyZero then 

begin 
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Dif.Re := 0; 
Dif.im := 0; 
end 
else 
Divide(Denom, Numer, Dif); { The difference is the } 
{ inverse of the fraction } 
end; { procedure ConstructDifference } 


function TestForRoot(X, Dif, Y, Tol : Extended) : boolean; 


{= These are the stopping criteria. Four different ones are 

{= provided. If you wish to change the active criteria, simply 

{= comment off the current criteria (including the appropriate OR) 
{= and remove the comment brackets from the criteria (including 
{= the appropriate OR) you wish to be active. 


{feeeeee eee ee eee ee ee ee ee ee eee eee ee ee eee 
begin 
TestForRoot := {-------- --------- } 
(ABS(Y) <= TNNearlyZero) {- Y=0 -} 
{- -} 
or {- -} 
+ 
(ABS(Dif) < ABS(X * Tol)) {- Relative change in X___ -} 
iE | 
{- } 
(Ol *) ie -} 
(* 2) { 1 
(* (ABS(Dif) < Tol) *) {- Absolute change in X___ -} 
(° i) if } 
("oF ‘ ( } 
(G 2) e a 
(* (ABS(Y) <= Tol) *) {- Absolute change in Y_ -} 


= The first criteria simply checks to see if the value of the 
= function is zero. You should probably always keep this criteria 
= active. 


= The second criteria checks the relative error in X. This criteria 
= evaluates the fractional change in X between interations. Note 

= that X has been multiplied throught the inequality to avoid divide 
= by zero errors. 


il 
Oe ee ne ee a ee | 


a ae a See Se 


The third criteria checks the absolute difference in X between 
{= iterations. 


tm | 


= The fourth criteria checks the absolute difference between 
the value of the function and zero. 


end; { procedure TestForRoot } 


begin { procedure FindOneRoot } 
Root := GuessRoot; 
Found := false; 
Iter := Q; 
EvaluatePoly(Degree, Poly, Root, yValue, yPrime, yDoublePrime); 
while (Iter < Maxiter) and not(Found) do 
begin 
Iter := Succ(lter); 
ConstructDifference(Degree, yValue, yPrime, yDoublePrime, Dif); 
Sub(Root, Dif, Root); 
EvaluatePoly(Degree, Poly, Root, yValue, yPrime, yDoublePrime); 
Found := TestForRoot(Modulus(Root), Modulus(Dif), Modulus(y Value), Tol); 
end; 
if not(Found) then Error := 1; { Iterations execeeded Maxlter } 
end; { procedure FindOneRoot } 


procedure ReducePoly(var Degree : integer; 
var Poly : TNCompVector; 
Root : TNcomplex); 


{= Input: Degree, Poly, Root 
{= Output: Degree, Poly 


{= This procedure deflates the polynomial Poly by 
{= factoring out the Root. Degree is reduced by one. 


var 
Term : integer; 
NewPoly : TNCompVector; 
Dummy : TNcomplex; 


begin 
NewPoly[Degree - 1] := Poly[Degree]; 
for Term := Degree - 1 downto 1 do 
begin 
Mult(NewPoly[Term], Root, Dummy); 
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Add(Dummy, Poly[Term], NewPoly[Term - 1]); 
end; 
Degree := Pred(Degree); 
Poly := NewPoly; 
end; { procedure ReducePoly } 


begin { procedure Laguerre } 
InitAndTest(Degree, Poly, Tol, Maxlter, InitGuess, NumRoots, Roots, 
yRoots, Iter, GuessRoot, InitDegree, InitPoly, Error); 
while (Degree > 0) and (Error = 0) do 
begin 
FindOneRoot(Degree, Poly, GuessRoot, Tol, Maxlter, 
Roots[NumRoots + 1], yRoots[NumRoots + 1], 
lter[NumRoots + 1], Error); 
if Error = 0 then 


{= The next statement refines the approximate root by = 
{= plugging it into the original polynomial. This = 
{= eliminates a lot of the round-off error = 
{= accumulated through many iterations = 


FindOneRoot(InitDegree, InitPoly, Roots[NumRoots + 1], 
Tol, Maxliter, Roots[NumRoots + 1], 
yRoots[NumRoots + 1], Addlter, Error); 
lier[NumRoots + 1] := lter[NumRoots + 1] + Addlter; 
NumRoots := Succ(NumRoots); 
ReducePoly(Degree, Poly, Roots[NumRoots]); { Reduce polynomial } 
end; 
GuessRoot := Roots[NumRoots]; 
end; 
end; { procedure Laguerre } 


begin 
end. { RootsOfEquat } 
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unit Message (6000); 


eee | ee eee ee ee SS SS Se eS ee 
= This unit provides the several messages which inform the 
{= user with some warnings and the help informations. 
{fesse 555225552 5555525555255 = 


{$T APPLFFTD} { Set the application type and creator } 
{$S+} 

{$l-} { Turn off I/O error checking } 

{$U RootsFinder} 

{$U SpecVar} 

{$U TurboGraph} 


interface 
uses 
MemTypes, QuickDraw, OSIntf, Toollntf,PackIntf, 
PasPrinter,SANE,MacPrint, RootsFinder, SpecVar, 
{$S Second Segment} 
TurboGraph; 


procedure SetUpWindows; 
procedure MakelnfoScreen; 
procedure DoAbout; 

procedure AlertBox1; 

procedure AlertBox2; 

procedure AlertBox3; 

procedure InfoGetEQParameter; 
procedure InfoGetCoeff; 
procedure InfoPlotOneParameter; 
procedure InfoPlotfwoParameter; 
procedure InfoPrint; 


implementation 


procedure SetUpWindows; 


begin 
{ The initial information screen } 
DefineWindow(InfoScrnWind, 100, 100, 362, 200, dBoxProc); 


{ The real transform window } 
DefineWindow(RootLocusWind, XMinGlb+5, YMinGlb+43, 
documentProc); 
DefineHeader(RootLocusWind, ‘Root Locus’); 


a1 


XMaxGlb-5, YMaxGlb-5, 


DefineWindow(RegendBox,100, 100,362,150, dBoxProc); 


{ The about box window } 

DefineWindow(AboutBoxWind, 125, 60, 387, 272, dBoxProc); 
DefineWindow(AlertBox,50, 100,360,170, dBoxProc); 
Define Window(HelpWind,20,60,480,330,dBoxProc); 


end; { SetUpWindows } 


procedure MakelnfoScreen; 


var 
SavePort : GrafPtr; 
InLeft : integer; { Window offset from left } 
InTop _: integer; { Window offset from top } 


begin 

GetPort(SavePort); { save the current grafport } 
selectWind(InfoScrnWind, FALSE); 

with Window([!nfoScrnWind].P’.portRect do 
begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 

end; 

MoveWindow(Window(InfoScrnWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(InfoScrnWind, TRUE); 
TextFont(SystemFont); 

TextSize(1 4); 

TextStyle([]); 

MoveTo(20, 20): 

DrawString(‘Welcome To’); 

TextSize(16); 

MoveTQ(70, 50); 

DrawString((MAC RootLocus’); 

TextSize(12); 

TextStyle([]); 

McveTo(90, 80); 

DrawString(‘Have a fun 1!'); 

repeat until Button; 

Hide Window( Window([InfoScrnWing].P); 
GetPort(SavePort); { save the current grafport } 

end; { MakelnfoScreen } 
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procedure DoAbout:; 


var 
SavePort : GrafPtr; 
InLeft : integer; { Window offset from left } 
InTop _: integer; { Window offset from top } 


begin 
GetPort(SavePort); { save the current grafport } 
SelectWind(AboutBoxWind, FALSE); 
with Window[AboutBoxWind].P’.portRect do 
begin 
InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 
MoveWindow(Window/[AboutBoxWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(AboutBoxWind, TRUE); 
TextFont(SystemFont); 
TextSize(18); 
TextStyle((]); 
MoveTo(60, 40); 
DrawString(‘(MAC RootLocus’); 
TextStyle({]); 
TextSize(12); 
MoveTo(85, 60); 
DrawString(‘version 1.00’); 
MoveTo(35, 90); 
Drawstring(‘Computer Aided Desgin Tool’): 
MoveTo(35, 110); 
DrawString('For RootLocus Analysis’): 
MoveTo(97, 160); 
DrawString(‘written by’); 
TextSize(13); 
MoveTo(85, 185); 
DrawString((KO SUNG HOON’); 
TextFont(SystemFont): 
repeat until Button; 
HideWindow(Window[AboutBoxWindg].P); 
SetPort(SavePort); 
end; { DoAbout } 


procedure AlertBox1; 


var 
SavePort : GrafPtr; 


InLeft : integer; { Window offset from left } 
InTop _: integer; { Window offset from top } 
begin 


GetPort(SavePort); { save the current grafport } 
SelectWind(3, FALSE); 


with Window([3].P*.portRect do 
begin 


inLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 


InTop := (YMaxGlib - (Bottom-Top)) div 2; 
end; 


MoveWindow(Window{[3]}.P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(3, TRUE); 


TextFont(SystemFont); 

TextSize(12); 

TextStyle((]); 

MoveTo(30, 30); 

DrawString(‘There is no Initial Degree.’); 
repeat until Button; 


HideWindow(Window{[3].P); 
end; 


procedure AlertBox2; 


[2s 55555555555 5555555555555555555==5) 
{= Bring up a window with a description of warning message = 
{2-2 eee wee eee ee eee eee sees sees ese=5} 
var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop __: integer; { Window offset from top } 
begin 


GetPort(SavePort); { save the current grafport } 
SelectWind(5, FALSE); 


with Window[5}.P’.portRect do 
begin 


InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 


InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 


MoveWindow(Window([5].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(5, TRUE); 


TextFont(SystemFont); 
TextSize(12); 
TextStyle([]); 
MoveTo(30, 30); 
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DrawString(‘There is no Initial Degree or’); 
MoveTo(50,50); 
DrawString(‘Characteristic Equation Coefficients.'); 
repeat until Button; 
Hide Window(Window([5].P); 

end; 


procedure AlertBox3; 


var 
SavePort : GrafPtr; 
InLeft : integer; { Window offset from left } 
InTop  : integer; { Window offset from top } 
begin 
GetPort(SavePort); { save the current grafport } 
SelectWind(5, FALSE); 
with Window[5].P’.portRect do 
begin 
InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 
MoveWindow(Window[5].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(5, TRUE); 
TextFont(SystemFont); 
TextSize(12); 
TextStyle([]); 
MoveTo(30, 30); 
DrawString(’You have a wrong Input.'); 
MoveTo(50,50); 
DrawString(‘'Check the Help Menu and try again.'); 
repeat until Button; 
HideWindow(Window{[5].P); 
end; 


procedure InfoGetEQParameter; 


var 
SavePort : GrafPtr; 
InLeft : integer; { Window offset from left } 
InTop _: integer; { Window offset from top } 


begin 
GetPort(SavePort); { save the current grafport } 
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SelectWind(HelpWind, FALSE); 
with Window[HelpWing].P’.portRect do 
begin 
InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 
MoveWindow(Window[HelpWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(HelpWind, TRUE); 
TextFont(SystemFont); 
TextSize(14); 
TextStyie({}); 
MoveTo(150, 15); 
DrawString’EQ PARAMETER); 
TextStyle((]); 
TextSize(12); 
MoveTo(20, 40); 
DrawString(’The EQ Parameter stands for equation parameter. It will allow you’); 
MoveTo(10, 55); 
DrawString(‘to input the degree of polynomial and some equation parameters. '); 
MoveTo(10,70); 
DrawString(’The degree of the polynomial should be 1 up to 10. Then there are ’); 
MoveTo(10, 85); 
DrawString(‘the default values for the other parameters. These values avoid '); 
MoveTo(10, 100); 
DrawString(’the convergence error for almost all polynomials. But if '); 
MoveTo(10, 115); 
DrawString(‘convergence error messages appears on the screen, you can change ’); 
MoveTo(10, 130); 
DrawString(‘these parameter values. These parameters must Satisfy the '); 
MoveTo(10, 145); 
DrawString(following conditions:’); 
MoveTo(30, 160); 
DrawString('(1) Initial guess >= 0 '); 
MoveTo(30, 175); 
DrawString('(2) Maximum Iteration >= 0 ’); 
MoveTo(30, 190); 
DrawString('(3) Tolerance > 0‘); 
MoveTo(40, 260); 
DrawString(” CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU ”); 
repeat until Button; 
Hide Window(Window[HelpWing}.P}; 
SetPort(SavePort); 
end; 


procedure InfoGetCoeff; 


96 


SavePort : GrafPtr; 
InLeft =: integer; { Window offset from left } 
InTop =: integer; { Window offset from top } 


begin 
GetPort(SavePort); { save the current grafport } 
SelectWind(HelpWind, FALSE); 
with Window[HelpWind].P’.portRect do 
begin 
InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 
MoveWindow(Window([HeloWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(HelpWind, TRUE); 
TextFont(SystemFont); 
TextSize(14); 
TextStyle([]); 
MoveTo(170, 15); 
DrawString('(GET COEFF’); 
TextStyle([]}); 
TextSize(12); 
MoveTo(20, 40); 
DrawString(‘The Get Coeff stands for Get Coefficients. It will allow you to’); 
MoveTo(10, 55); 
DrawString(‘input the algebraic expression for the coefficients of the  '); 
MoveTo(10,70); 
DrawString(‘characteristic equation. It may have up to two undetermined '); 
MoveTo(10, 85); 
DrawString(‘parameters(A and B). In case of the one parameter root locus '); 
MoveTo(10, 100); 
DrawString(’method, you use only one undetermined paramerter(A). The routine’); 
MoveTo(10, 115); 
DrawString(‘uses standard algebraic, or infix, notation with parenthesis’); 
MoveTo(10, 130); 
DrawString(‘allowed. Operater can include +, -, *, /, and “(expomentiation).’); 
MoveTo(10, 145); 
DrawString(‘The unary minus sign is allowed.'); 
MoveTo(10, 160); 
DrawString(‘If you choose the Get Coeff command without the degree of the’); 
MoveTo(10, 175); 
DrawString(‘polynomial,the message appears on the screen. It tells you that’); 
MoveTo(10, 190); 
DrawString(‘the degree of the polynomial has not been entered yet.’); 
MoveTo(40, 260); 
DrawString("" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU "); 


wf 


repeat until Button; 
Hide Window(Window[HelpWind].P); 
SetPort(SavePort); 

end; 


procedure InfoPlotOneParameter; 


var 
SavePort : GrafPtr; 
InLeft : integer; { Window offset from left } 
InTop __: integer; { Window offset from top } 


begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(HelpWind, FALSE); 

with Window[HelpWind].P’.portRect do 
begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[HelpWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(HeipWind, TRUE); 

TextFont(SystemFont); 

TextSize(14); 

TextStyle({]); 

MoveTo(145, 15); 

DrawString((ONE PARAMETER’); 

TextStyle([]); 

TextSize(12); 

MoveTo(20, 35); 

DrawString(‘One parameter Command will allow you to input the plot data of’); 
MoveTo(10, 50); 

DrawString(’the one parameter root locus method. The default values are shown’); 
MoveTo(10,65); 
DrawString(‘in the dialog box. They can be changed as desired.'); 

MoveTo(290, 80); 

DrawString('-. Enter the minimum and maximum gain values into the Min Gain’); 
MoveTo(10, 95); 

DrawString(‘and the Max Gain insertion box.’); 

MoveTo(20, 110); 

DrawString(’-. Choose the type of interval. There are two types of interval:’); 
MoveTo(10, 125); 

DrawString(’Linear and Logarithmic interval. When you click the radio button,’); 
MoveTo(10, 140); 

DrawString(‘the desired type of interval is choosen.’); 
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MoveTo(20, 155); 
DrawString('-.Choose the type of scale for the axis. There are two types of ’); 
MoveTo(10, 170); 
DrawString(’scale: Auto and Manual scale. When you click the radio button,’); 
MoveTo(10, 185); 
DrawString('the desired type of scale is choosen.'); 
MoveTo(20, 200); 
DrawString(’-.Enter some number into Points to Plot inserton box in order to '); 
MoveTo(10, 215); 
DrawString('‘decide the plot resolution.’); 
MoveTo(40, 260); 
DrawString("" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU "; 
repeat until Button; 
HideWindow(Window[HelpWing].P); 
SetPort(SavePort); 
end; 


procedure InfoPlotTtwoParameter; 


var 
SavePort : GrafPtr; 
InLeft : integer; { Window offset from left } 
InTop _: integer; { Window offset from top } 


begin 
GetPort(SavePort); { save the current grafport } 
SelectWind(HelpWind, FALSE); 
with Window[HelpWind].P*.portRect do 
begin 
InLeft := (XMaxGlib - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 
MoveWindow(Window[HelpWing].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(HelpWind, TRUE); 
TextFont(SystemFont); 
TextSize(14); 
TextStyle((]); 
MoveTo(145, 15); 
DrawString('TWO PARAMETER’); 
TextStyle([{]); 
TextSize(12); 
MoveTo(20, 35); 
DrawString(‘Two parameter command will allow you to input the plot data of); 


MoveTo(10, 50); 
DrawString(‘the two parameter root locus method. The default values are shown), 
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MoveTo(10,65); 

DrawString(‘in the dialog box. They can be changed as desired.’); 

MoveTo(20, 80); 

DrawString(’The items in this dialog are similar to those in the one parameter’); 
MoveTo(10, 95); 

DrawString(‘plot data dialog box, but several items are different. ’); 

MoveTo(20, 110); 

DrawString(‘There exists one more undetemined parameter B to be inserted. The’); 
MoveTo(10, 125); 

DrawString('How many loci item lets you decide how many loci are to be drawn '); 
MoveToac10, 140); 

DrawSirina(‘for each parameter. It will be 1 up to 10 for each parameter.’); 
MoveTo(10, 155); 

DrawString(‘There is no auto scale for axes. Only the manual scale is available’); 
MoveTo(10, 170); 

DrawString(The last item is the marking and justification in order to draw’); 
MoveTo(10, 185); 

DrawString(‘the selected A and B values on the plot. Two buttons are chosen’); 
MoveTo(10, 200); 

DrawString(‘each time for each parameter; one for position, the other justifi- ‘); 
MoveTo(10, 215); 

DrawString(‘cation to draw."); 

MoveTo(40, 260); 
DrawSiring("™ CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU "); 

reoeat until Button; 

Hide Window(Window[HelpWind].P); 

setPort(SavePort); 

end; 


procedure InfoPrint; 


var 
SavePort : GrafPtr; 
InLeft : integer; { Window offset from left } 
InTop  : integer; { Window offset from top } 


begin 
GetPort(SavePort); { save the current grafport } 
SelectWind(HelpWind, FALSE); 
with Window[HelpWind]}.P’.portRect do 
begin 
InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 
MoveWindow(Window[HelpWing].P, InLeft, InTop, TRUE); { Center the window } 
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SetVisibility(HelpWind, TRUE); 

TextFont(SystemFont); 

TextSize(14); 

TextStyle({]); 

MoveTo(160, 15); 

DrawString(‘PRINT OUT'); 

TextStyle({]); 

TextSize(12); 

MoveTo(20, 35); 

DrawString(‘In order to print out your work, there are a few ways available.’); 
MoveTo(20, 50); 

DrawString(’-.Select the print command in the File menu. This allows the '); 
MoveTo(10,65); 

DrawString(‘user to get a hard copy of any plot displayed.’); 

MoveTo(20, 80); 

DrawString('-. Hold the command key and then type the number 4 to print '); 
MoveTo(10, 95); 

DrawString(‘the contents of the active window immediately.’); 

MoveTo(20, 110); 

DrawString('-. Press the command and shift key and type the number 3 to '); 
MoveTo(10, 125); 

DrawString(‘create a MacPaint document.’); 

MoveTo(40, 260); 
DrawString(" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU "); 
repeat until Button; 

HideWindow(Window[HelpWing].P); 

SetPort(SavePort); 

end; 


begin 
end.{unit message} 
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unit MyDialog (7000); 


{= This unit supports four dialog boxes for input data. =} 
Bo 5 5 esa 
{$B+} { Set the bundle bit } 

{$R MacRootLocus.Rsrc} { Identify resource file for menu and icon info } 
{$T APPLFFTD} { Set the application type and creator } 

{$S+} { generate segmented code} 

{$l-} { Turn off 1/O error checking } 


{SU SpecVar} 
{$U RootsFinder} 
{$U MakeRoot} 
{$U TurboGraph} 
{$U Message} 


interface 
uses 
MemTypes, QuickDraw, OSIntf, Toollntf,Packintf, 
PasPrinter,SANE,MacPrint, RootsFinder, SpecVar, 
{$S Second Segment} 
TurboGraph, 
Message, 
{$S Third Segment} 
MakeRoot; 


procedure MakeRegend; 
procedure GetEQParameter; 
procedure GetCoeff; 

procedure PlotOneParameter; 
procedure PLotTwoParameter; 


implementation 


procedure MakeRegend; 


— This procedure provides the informaton box with =} 
= the capabilety of the text editor =} 
{fe s==S=====5=25=25 2525255252222 2222222222 =} 
var 


SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 
InTop __: integer; { Window offset from top } 
begin 

GetPort(SavePort); { save the current grafport } 
SelectWind(3, FALSE); 
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Textsize(2); 

with Window(3].P*’.portRect do 

begin 
InLeft : 


InTop := Round((YMaxGlb - (Bottom-Top)) /1.05); 
end; 


Round((XMaxGlb - (Right-Left)) / 1.38); { Calculate window offsets } 


MoveWindow(Window[3].P, InLeft, InTop, TRUE); { Center the window } 


SetVisibility(3, TRUE); 
txRect := thePort’.portRect; 
textH := TENew(txRect, txRect); 
TEldle(textH); 
TextinputEnabled := True; 
end; 


procedure GetEQParameter; 


{ This procedure provides the dialog box for the 
{= equation parameter of the characteristic equation. 
{e=2=2=25222255555222552522222222222=5 
const 

CancelBtn = 1; 

OKBtn = 2; 

InitDegreeText = 3; 

Inita@uessText = 4; 

MaxlterText = 2) 

ToleranceText = 6; 

var 


tempdegree,templter : Longint; 
tempGuess : TNcomplex; 
tempTolerance : Extended; 

BoxPlot : Boolean; 
begin 

theDialog := GetNewDialog(256, nil, pointer(-1)); 
GetDitem(theDialog, 3, theType, h3, r); 
GetDitem(theDialog, 4, theType, h4, r); 
GetDitem(theDialog, 5, theType, h5, r); 
GetDlitem(theDialog, 6, theType, h6, r); 
Done := False; 
tempGuess.Re := 1.0; 
tempGuess.|m := 0.0; 

templter := 100; 
tempTolerance := 1E-7; 
BoxPlot := False; 
repeat 

ModalDialog(nil, itemHit); 

case itemHit of 

CancelBtn : done := True; 


103 


OKBtn : begin 
InitDegree := tempdegree; 
InitGuess.Re := tempGuess.Re; 
InitGuess.Im := tempGuess.|m; 
Maxiter := templter; me 
Tolerance := tempTolerance; 
done ‘= True; 
end; 
InitDegreeText : begin 
GetiText(h3, s); 
StringTaNum(s, tempdegree); 
InitDegreeStatus := True; 
if (tempdegree <1 ) or (tempdegree >10) then 
begin 
sysbeep(10); 
end; 
end; 
InitauessText : begin 
GetlText(h4, s); 
tempGuess.Re := Str2Num(s); 
if (tempGuess.Re <0) then 
begin 
sysbeep(10); 
end; 
end; 
MaxlterText : begin 
GetlText(h5, s); 
StringToNum(s, templter); 
if (templter <0 ) then 
begin 
sysbeep(10); 
end; 
end; 
ToleranceText : begin 
GetlText(h6, s); 
tempTolerance := Str2Num(s); 
if (tempTolerance <= 0) then 
begin 
sysbeep(10); 
end; 
end; 
end; 
until done; 
DisposDialog(the Dialog); 
end; 


procedure GetCoeff; 
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{ This procedure provides the dialog box for the =} 
= coefficients of the characteristic equation. =} 
=eSS55555555555555 5555525222225} 
type char_set = set of char; 
const 

CancelBtn =1; 

OKBtn = 2; 

SoText = 3; 

S1Text = 4; 

S2Text = 5; 

S3Text = 6; 

S4Text = 7; 

S5Text = 8; 

Sé6Text =r 

S7Text = 10; 

S8Text = 11; 

S9oText = 12; 

$10Text = 13; 

var 

Idno,tempitemhit,sssindex,ssslength 

: Integer; 
h : Array[1..13] of Handle; 


Sset_valid :char_set; 

op_set : char_set; 

sss: str255; 

sssc: char; 
noerrorcheck,leftparen: boolean; 


Procedure skipblank ; 
begin 
while ((sss[sssindex] = '’) and (sssindex <= ssslength)) do 
sssindex := sssindex + 1; 
end; 


Procedure Checknext ( setofnext1,setofnext2:char_set); 
var opchar: char; 
begin 
opchar := SSSsC; 
sssindex := sssindex +1; 
SSSC := sss[sssindex]; 
if (sssc =’ ’) then begin 
skipblank; 
SSSC := sss[sssindex]; 
if sssindex <= ssslength then begin 
if not (sssc in setofnext1) then 
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noerrorcheck:= False; 
end else 
if (opchar in (op_set + ['('])) then 
noerrorcheck := false 
end else 
if not (sssc in setofnext2) then 
noerrorcheck:= false 
end; 


Procedure CheckExpression(matchparen:boolean); 
begin 
lefiparen:= true; 
while ((noerrorcheck and leftparen) and (sssindex < ssslength)) do 
case (sssc) of 
’0'..'9': Checknext(op_set+['.’,’)'],op_set+['0"..'9',".",')']); 
‘a',b': Checknext(op_set+ [')'],op_set+[‘a’,'b’,’)']); 
rer rary t-' 4's Checknext(['0'’..'9','a','b','(‘},['0'..'9','a’,'b', (J); 
"": Checknext(['0'..'9'],['0'..'9']); 
‘(begin 
Checknext(['0'..'9','a’,'b',’-','(‘],['0"..’9','a', ‘b's, (i 
CheckExpression (false); 
if sssc = ')' then 
if sssindex < ssslength then begin 
SSSC ‘= $$s{sssindex]; 
if sss[sssindex + 1] =’’ then skipblank; 
if sssindex < ssslength then 
Checknext(op_set+[')'],op_set+]')']) 
else if not matchparen then noerrorcheck := false 
end else begin 
if not matchparen then noerrorcheck := false 
end 
else noerrorcheck := false; 
end; 
')':begin 
if matchparen then noerrorcheck:= false; 
leftparen := false; 
end; 
‘begin 
skipblank ; 
if sssindex = ssslength then begin 
if not (sss[sssindex] in ['0'..'9','a’,'b‘]) then 
noerrorcheck:= false 
end else if sssindex < ssslength then sssc := sss[sssindex]; 
end 
end; 


if not noerrorcheck then sysbeep (30); 
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end; 
{begin 
if not InitDegreeStatus then 
begin 
SysBeep(30); 
AlertBox1; 
end 
else} 
begin 
set_valid := 


mL Omen a, D, te) sy eo eae 
Oeeset i= ['+',/',-'7°',4'); 
if not InitDegreeStatus then 
begin 
SysBeep(30); 
AlertBox1; 
end 
else 
begin 
Idno := 300 + InitDegree; 
theDialog := GetNewDialog(Idno, nil, Pointer (-1)); 
for n := 3 to InitDegree + 3 do 


GetDltem(theDialog, n, theType, h[n], 1); 
done := False; 


ModalDialog(nil, itemHit); 


repeat 


tempitemhit := itemHit; 
case itemHit of 


CancelBtn : done := True; 
OKBtn : done := True; 
3.12 : begin 


while tempitemhit = itemHit do 
ModalDialog(nil, tempitemhit); 
GetlText (h[itemHit],sss); 
noerrorcheck := true; 
sssindex:= 1; 
ssslength := Length (sss); 
SSsc ‘= sss[1]; 
if not (sss[1] in ['0'..'9','a’,'b’,-','(‘J) then 
noerrorcheck := False; 
CheckExpression ( True); 
GetlT ext(h[itemHit],InfixArray[InitDegree + 3 - itemHit)); 
GetCoeffStatus := True; 
itemhit := tempitemhit; 
if itemhit = OKBtn then done := true 
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end; 
end; 
until done; 
DisposDiaiog(theDialog); 
end; 
end; 
{end;} 


procedure PlotOneParameter; 


{ 

{= This procedure provides the dialog box for the =} 
= plot data of the one parameter root locus method. =} 
{ 


CancelBtn = 
OKBtn = 2; 
LinearBtn = 3; 

LogarithmicBtn = 4; 
AutoBtn = 5; 

ManualBtn = 6; 
AMinGainText =7; 
AMaxGainText =8 
XMinT ext =9 
XMaxText = | 
YMinText = | 
Y¥MaxtText = | 
PlotPointsText = 


var 
saveSoundVol : Integer; 
radButton : array [3..6] of ControlHandle; 
h : Array[7..13] of Handie; 
tempAMinGain : Extended; 
tempAMaxGain_: Extended; 
tempXMin : Extended; 
tempXMax : Extended; 
tempYMin : Extended; 
temp Max ; Extended; 
tempPoints : Longlnt; 
OkPiot ‘boolean; 
begin 
if not InitDegreeStatus or not GetCoeffStatus then 
begin 
SysBeep(30); 
AlertBox2; 
end 
else 
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begin 
GetSound\ol (saveSoundVol); 
SetSoundVol (1); 
theDialog := GetNewDialog(257, nil, Pointer (-1)); 
for n := LinearBtn to ManualBtn do 
GetDltem(theDialog,n,theType,Handle(radButton{n)),r); 


SetCtlValue (radButton{LinearBtn],1); 
linear := True; 


setCtlValue (radButton{AutoBtn],1); 
AutoScale := True; 
for n := 7 to 13 do 
GetDlitem(theDialog, n, theType, h{n}, r); 
tempAMinGain := Str2Num('0.1’); 
tempAMaxGain := Str2Num('10000’'); 
tempXMin := Str2Num(’-10'); 
tempXMax := Str2Num('5'); 
tempYMin := Str2Num('-10'); 
tempYMax := Str2Num('10’); 
StringToNum('50', tempPoints); 
done := False; 
OkPlot := False; 
repeat 
ModalDialog(nil,itemHit); 
case itemHit of 


CancelBtn : done := True; 
OKBtn : begin 
done := True; 
OkPlot := True; 
end; 
LinearBtn : begin 


SetSoundVol (1); 


for n := LinearBtn to LogarithmicBtn do 


SetCtlValue(radButton[n},Ord(n = itemHit)); 
linear := True; 
end; 


LogarithmicBtn : begin 
SetSoundVol (1); 
for n := LinearBtn to LogarithmicBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
linear := False; 
end; 
AutoBtn : begin 
SetSoundVol (1); 
for n := AutoBtn to ManualBtn do 


SetCtlValue(radButton[(n],Ord(n = itemHit)); 
AutoScale := True; 
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end; 


ManualBtn : begin 
SetSoundVol (1); 
for n := AutoBtn to ManualBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
AutoScale := False; 
end; 
AMinGainText  : begin 
GetlText(h[itemHit] , s); 
tempAMinGain := Str2Num(s); 
if (tempAMinGain < 0) or (tempAMinGain > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
AMaxGainText  : begin 
GetiText(h{itemHit] , s); 
tempAMaxGain := Str2Num(s); 
if (tempAMaxGain <0) or (tempAMaxGain > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
XMinText : begin 
GetiText(hf{itemHit] , s); 
tempXMin := Str2Num(s); 
if (tempXMin <-1e7 ) or (tempXMin > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
XMaxText : begin 
Get!Text(h[itemHit] , s); 
tempXMax := Str2Num(s); 
if (tempXMax <-1e7 ) or (tempXMax > 1e7) then 
begin 
sysbe2p(10); 
end; 
end; 
YMinText : begin 
GetlText(h[itemHit] , s); 
tempYMin := Str2Num(s); 
if (tempYMin <-1e7 ) or (tempYMin > 1e7) then 
begin 
sysbeep(10); 
end; 
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end; 
YMaxText : begin 
GetlText(h{itemHit] , s); 
tempYMax := Str2Num(s); 
if (tempYMax <-1e7 ) or (tempYMax > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
PlotPointsText : begin 
GetlText(h[itemHit] , s); 
StringToNum(s, tempPoints); 
if (tempPoints <= 0) or (tempPoints > 150) then 
begin 
sysbeep(10); 
end; 
end; 


end;{case end} 

until done; 
DisposDialog(theDialog); 

if OkPlot then 

begin 

AMinGain := tempAMinGain; 
AMaxGain := tempAMaxGain; 


XMn = tempXMin; 
XMx := tempXMax; 
YMn = tempYMin; 
YMx = tempYMax; 
Points := tempPoints; 
HideCursor; 
Rewrite(OutFile,,RootLocus1 .data’); 

GetRoott; 
Close(OutFile); 
ShowCursor; 

MakeRegend; 

end;{if} 

end;{case} 


end:{PlotOneParameter} 


procedure PLotTwoParameter; 


Ss This procedure provides the dialog box for the = 
& plot data of the two parameter root locus method. = 
{o=255555525252255525555555555555555) 
const 

CancelBtn = 1; 


AA 


OKBtn = 2; 


LinearBtn ; 
LogarithmicBtn = 4; 
AMarkStBtn = 55 
AMarkEadBtn = 6; 
AJustifyRBtn = 7; 
AJustifyLBtn = 8; 
BMarkStBtn = 9; 
BMarkEdBtn = 10; 
BJustifyRBin = 11; 
BJustifyl Btn = 12; 
AMinGainText = 13; 


AMaxGainText = 14; 
BMinGainText = 15; 
BManGainText = 16; 


LociText = 17; 
PlotPointsText = 18; 
XMinText = 19; 
XMaxText = 20; 
YMinText = 21; 
YMaxText = 22; 
var 


saveSoundVol,i : Integer; 
radButton : array [3..12] of ControlHandle; 
h : Array[13..22] of Handle; 
tempAMinGain : Extended; 
tempAMaxGain _: Extended; 
tempBMinGain : Extended; 
tempBMaxGain_ : Extended; 
tempXMin : Extended; 
tempXMax : Extended; 
tempYMin : Extended; 
tempYMax : Extended; 
tempStep, tempPoints : Longint; 
OkPlot : Boolean; 
begin 
if not initDegreeStatus or not GetCoeffStatus then 
begin 
SysBeep(30); 
AlertBox2; 
end 
else 
begin 
GetSoundVol (saveSoundVol); 
SetSoundVol (1); 
theDialog := GetNewDialog(258, nil, Pointer (-1)); 


2 





for n := LinearBtn to BJustifyLBtn do 
GetDlitem(theDialog,n,the Type, Handle(radButton[n)}),r); 
setCtlValue (radButton[LinearBtn],1); 
linear := True; 
setCtlValue (radButton[AMarkStBtn],1); 
AMarkStatus := False; 
SetCtlValue (radButton[AJustifyRBtn],1); 
ARJustification := True; 
SetCtlValue (radButton[BMarkStBtn],1); 
BMarkStatus := False; 
setCtlValue (radButton[BJustifyRBtn],1); 
BRuJustification := True; 
for n := 13 to 22 do 
GetDlitem(theDialog, n, theType, h[n], r); 
tempAMinGain := Str2Num('0.1'); 
tempAMaxGain := Str2Num('10000'); 
tempBMinGain := Str2Num('0.1'); 
tempBMaxGain := Str2Num('10000'); 
tempXMin := Str2Num('-10'); 
tempXMax := Str2Num('5'); 
tempYMin := Str2Num('-10'); 
tempYMax := Str2Num('10'); 
StringToNum('4', tempStep); 
StringToNum('50', tempPoints); 
done := False; 
OkPlot := False; 
repeat 
ModalDialog(nil,itemHit); 
case itemHit of 
CancelBtn :done := True; 
OKBtn : begin 
done := True; 
OkPlot := True; 
end; 
LinearBtn : begin 
SetSoundVol (1); 
for n := LinearBtn to LogarithmicBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
linear := True; 
end; 


LogarithmicBtn : begin 
SetSoundVol (1); 
for n := LinearBtn to LogarithmicBtn do 
SetCtiValue(radButton[n],Ord(n = itemHit)); 
linear := False; 
end; 


is 


AMarkStBtn : begin 
SetSoundVol (1); 
for n := AMarkStBtn to AMarkEdBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
AMarkStatus := False; 
end; 


AMarkEcBtn : begin 
SetSoundVol (1); 
for n := AMarkStBtn to AMarkEcBtn do 
SetCiiValue(radButton[n],Ord(n = itemHit)); 
AMarkSitatus := True; 
end; 
AJustifyRBtn : begin 
SetSoundVol (1); 
for n := AJustifyRBtn to AJustifyLBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
ARJustification := True; 
end; 


AJustifyLBtn : begin 
SetSoundVol (1); 
for n := AJustifyRBtn to AJustifyLBtn do 
SeiCiValue(radButton[n],Ord(n = itemHit)); 
ARJustification := False; 
end; 
BMarkStBtn : begin 
SetSoundVol (1); 
for n := BMarkStBtn to BMarkEdBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
BMarkStatus := False; 
end; 


BMarkEaBtn : begin 
SetSoundVol (1); 
for n := BMarkStBitn to BMarkEdBtn do 
SetCilValue(radButton[n],Ord(n = itemHit)); 
BMarkStatus := True; 
end; 
BJustifyRBtn : begin 
SetSoundVol (1); 
for n := BJustifyRBtn to BJustifyLBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
BRJustification := True; 
end; 


BJustifyLBtn : begin 
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SetSoundVol (1); 
for n := BJustifyRBtn to BJustifyLBtn do 
SetCtiValue(radButton[n],Ord(n = itemHit)); 
BRuJustification := False; 
end; 
AMinGainText : begin 
GetiText(h[itemHit],s); 
tempAMinGain := Str2Num(s); 
if (tempAMinGain < 0) or (tempAMinGain > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
AMaxGainText : begin 
GetiText(h[itemHit],s); 
tempAMaxGain := Str2Num(s); 
if (tempAMaxGain < 0) or (tempAMaxGain > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
BMinGainText : begin 
GetlText(h[itemHit],s); 
tempBMinGain := Str2Num(s); 
if (tempBMinGain < 0) or (tempBMinGain > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 


BManGainText : begin 
GetlText(h[{itemHit],s); 
tempBMaxGain := Str2Num(s); 
if (tempBMaxGain < 0) or (tempBMaxGain > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
XMinText : begin 
GetiText(h[itemHit],s); 
tempXMin := Str2Num(s); 
if (tempXMin <-1e7 ) or (tempXMin > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
XMaxText : begin 
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GetlText(h[itemHit],s); 
tempXMax := Str2Num(s); 
if (tempXMax <-1e7 ) or (tempXMax > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
YMinText : begin 
GetlText(h[itemHit],s); 
tempYMin := Str2Num(s); 
if (fampYMin <-1e7 ) or (tempYMin > 1e7) then 
begin 
sysbeep(10); 
end; 
end; 
YMaxText : begin 
GetlText(h[itemHit],s); 
tempYMax := Str2Num(s); 
if (tempYMax <-1e7 ) or (tempYMax > 1e7) then 
begin 
sysbeep(10);: 
end; 
end; 
LociText : begin 
Get!Text(hf{itemHit],s); 
StringToNum(s, tempStep); 
if (tempStep <= 0) or (tempStep > 15) then 
begin 
sysbeep(10); 
end; 
end; 
PlotPointsText : begin 
GetiText(h[itemHit],s); 
StringToNum(s, tempPoints); 
if (tempPoints <= 0) or (tempPoints > 150) then 
begin 
sysbeep(10); 
end; 
end; 
end:{case end} 
until done; 
DisposDialog(theDialog); 
if OkPlot then 
begin 
AMinGain — := tempAMinGain; 
AMaxGain := tempAMaxGain; 
BMinGain  := tempBMinGain; 
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BMaxGain := tempBMaxGain; 


XMn <= tempXMin; 
XMx = tempXMax; 
YMn = tempYMin; 
YMx = tempYMax; 
Step = tempStep; 
Points := tempPoints; 
HideCursor; 


Rewrite(OutFile,'RootLocus2.data’); 

GetRoot2; 

Close(OutFile); 

ShowCursor; 
MakeRegend; 

end;{if} 

end;{case} 
end;{Plot'woParameter} 


begin 
end.{unit MyDialog} 


unit TurboGraph(11000); 


= Turbo Pascal Numerical Methods Toolbox 
{= Copyright (C) 1987 Borland International 
{ 


= This unit provides routines for displaying graphics. = 


{ZU RootsFinder} 
{SU SpecVar} 


interface 


uses 
MemTypes, QuickDraw, OSintf, Toollntf, Packintf,PasPrinter, SANE, MacPrint , 
RootsFinder, SpecVar; 


const 
MaxWorldsGlb = 10; { The maximum number of worlds that can be defined } 
MaxWindowsGlb = 10; { The maximum number of windows that can be defined } 


{MaxPiotGlb = 1500;}{ The maximum number of points in a plot array } 
X1 Offset = 40; { The upper left corner of the axis } 
Y1Offset = 110) 
X2Offset = 30; { The lower right corner of the axis } 
Y2Offset = o0R 
type 


WrkString = Str255; { A general string type } 


WorldType = record { Used to store world coordinates } 
X1, Y1, X2, Y2 : real; 
end; 


WindowType = record { Used to store window information } 
X1, Y1, K2, Y2 : integer; { The windows screen coordinates } 
Header : WrkString; { Header for a window } 
W : WindowRecord; { Mac window record } 
P : WindowPtr; { Mac window pointer } 
H : PicHandle; { A handle to a picture } 
end; 


Worlds = array[1..MaxWorldsGlb] of WorldType; { Holds world info } 
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Windows = array[1..MaxWindowsGlb] of WindowType;{ Holds window info } 
{ PlotArray = array[1..MaxPlotGlb, 1..2] of Extended;} 


var 
{ Coordinates of the currently active window } 
XMinGlb, YMinGlb, XMaxGlb, YMaxGlb : integer; 


{ Coordinates of the currently selected world } 
X1WidGib, X2WIdGib, Y1WidGlb, Y2WIdGib : real; 


{ World coordinate scaling factors } 
AxGlb, AyGlb, BxGlb, ByGlb : real; 


{ Coordinates of an axis if one is defined } 
X1RefGlb, X2RefGlb, Y1RefGlb, Y2RefGlb : integer; 


{ The maximum world and window number defined } 
MaxWorldGlb, MaxWindowGlb : integer; 


{ The currently selected world and window } 
WorldNdxGib, WindowNdxGlb : integer; 


{ Aspect ratio for a true circle } 
AspectFactor : real; 
AspectGlb — : real; 


{ Flags if an axis is defined } 
AxisGlb : boolean; 


{ Flags if hatching is turned on } 
HatchGlb : boolean; 


{ Flags if clipping should be performed } 
ClippingGlb : boolean; 


{ The currently selected line style } 
LineStyleGlb : integer; 


{ The current foreground color } 
ForeColorGlb : integer; 


{ Holds the worlds defined by the user } 
World : Worlds; 


{ Holds the windows defined by the user } 
Window : Windows; 


ES) 


procedure Error(ProcName : WrkString); 
{ Report that an error occurred in the procedure named "ProcName” } 


procedure SetForegroundColor(Color : integer); 
{ Set the foreground drawing color } 


procedure SetBackgroundColor(Color : integer); 
{ Set the background color } 


procedure DP({X, Y : integer); 
{ Piot a pixel at position (A, Y) } 


function PD(X, Y : integer) : boolean; 
{ Return true if the color of the pixel at (X, Y) matches ForeColorGlb } 


procedure DrawStraight(X1, X2, Y : integer); 
{ Draw a horizontal line from X1,Y to X2,Y } 


procedure InvertWindow(WinNum : integer); 
{ Invert the window referenced by WinNum } 


function RealToStr(R : real) : WrkString; 
{ Returns the string representation of the real R. } 


function IntToStr(l : integer) : WrkString; 
{ Returns the string representation of the integer I. } 


procedure SetLineStyle(Style : integer); 
{ Select the current line style } 


function GetLineStyle : integer; 
{ Returns the current line style } 


procedure DefineWorld(WorldNum : integer; X_1, Y_1, X_2, Y_2 : real); 
{ Defines a world coordinate system with a specific number } 


procedure SelectWorld(WorldNum : integer); 
{ Select the worid associated with WorldNum } 


procedure DefineWindow(WinNum, XLo, YLo, XHi, YHi, WindowType : integer); 
{ Defines a window with a specific window number and window type } 


procedure SetVisibility(WinNum : integer; Visible : boolean); 
{ Sets the visibility of a window to either TRUE or FALSE } 


procedure SelectWind(WinNum : integer; Visible : boolean); 
{ Selects a window as visible or invisible } 


4 
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procedure DefineHeader(WinNum : integer; Hdr : WrkString); 
{ Define a header for a window } 


procedure RemoveHeader(WinNum : integer); 
{ Clears a header from a window } 


procedure ReDefineWindow(WinNum, X1, Y1, X2, Y2, WindowType : integer); 
{ Redefines the coordinates of an existing window } 


function WindowxX(X : real) : integer; 
{ Converts an X world coordinate into a screen coordinate } 


function WindowY(Y : real) : integer; 
{ Converts a Y world coordinate into a screen coordinate } 


procedure ResetWorlds; 
{ Resets all worlds to the maximum dimensions of the screen } 


procedure InitGraphic; 
{ Initializes the graphics system } 


function Clip(var X1, Y1, X2, Y2 : integer) : boolean; 
{ Clips a line to the current window and returns TRUE if any part } 
{ of the line is still in the window after the clip operation } 


procedure DrawPoint(Xr, Yr : real); 
{ Draw a point in world coordinates } 


function PointDrawn(Xr, Yr : real) : boolean; 
{ Returns TRUE if the point at (Xr, Yr) is drawn } 


procedure DrawLine(X1, Y1, X2, Y2 : real); 
{ Draw a line in world coordinates } 


procedure DrawLineClipped (X1, Y1, X2, Y2 : integer); 
{ Draw a line clipped in screen coordinates } 


procedure DrawCrossDiag(X, Y, Scale : integer); 
{ Draw across at screen coordinate (X, Y) with a scaling factor } 


procedure DrawWye(X, Y, Scale : integer); 
{ Draw a Y at screen coordinate (X, Y) with a scaling factor } 


procedure DrawDiamond(X, Y, Scale : integer); 
{ Draw a cross at screen coordinate (X, Y) with a scaling factor } 


procedure DrawCircleDirect(Xr, Yr, R : integer); 
{ Draw a circle at screen coordinate (Xr, Yr) with a radius R } 


procedure DrawCircle(X_R, Y_R, Xradius : real); 
{ Draw a circle at world coordinate (X_R, Y_R) with a radius R } 


procedure DrawCross(X1, Y1, Scale : integer); 
{ Draw a cross at screen coordinate (X, Y) with a scaling factor } 


procedure DrawStar(™, Y, Scale : integer); 
{ Draw a Star at screen coordinate (X, Y) with a scaling factor } 


procedure DrawSquareC(X1, Y1, X2, Y2 : integer; Fill : boolean); 
{ Draw a Square in screen Coordinates with optional filling } 


procedure DrawSquare(X1, Y1, X2, Y2 : real; Fill : boolean); 
{ Draw a square in world coordinates with optional filling } 


procedure DrawAscii(X, Y : integer; Size, CharByte : byte); 
{ Draw a character with ASCII code CharByte } 


procedure DrawText(X, Y, Scale : integer; Txt : Str255); 
{ Draw a string at screen coordinate (X, Y) with a scaling factor } 


procedure DrawTextW(X, Y : real; Scale : integer; Txt : WrkString); 
{ Draw a String at world coordinate (X, Y) with a scaling factor } 


procedure TextStyle(Face : Style); 
{ Face = (bold, italic, underline, outline, shadow, condense, extend) } 


procedure HardCopy(TopWin : boolean); 
{ Do a screen dump of the currently selected window (TopWin = TRUE) } 
{ or the whole screen (TopWin = FALSE) } 


procedure OpenPic(WinNum : integer; ShowPic : boolean); 
{ Open a picture for a specific window and only show the drawing } 
{ if SnowPic is set to TRUE } 


procedure ODrawPic(WinNum : integer); 
{ Draw the picture associated with a window } 


procedure ErasePic(WinNum : integer); 
{ Erase the picture associated with a window } 


procedure ClearWindow(WinNum : integer); 
{ Clear the content portion of a window } 


function Wherex : integer; 
{ Returns the X cursor position } 


function WhereyY : integer; 
{ Returns the Y cursor position } 


procedure SetWindow(X1, Y1, X2, Y2 : integer); 
{ Creates an invisible window inset from the currently selected one } 


procedure FindWorld(| : integer; A : PlotArray; NPoints : integer); 
{ Finds a world to fit a polygon defined in A } 


procedure FindWorldi(I : integer; XMn,YMn,XMx,YMx:ExTended ); 
{ Finds a world to fit a polygon defined in A } 


procedure DrawAxis(Footer1, Footer2 : WrkString; Arrows : boolean); 
{ Draws an axis with Footers and optional arrows on the axis } 


procedure ResetAxis; 
{ Sets AxisGlb to TRUE } 


procedure DrawPolygon(A : PlotArray; First, NPoints, Line, Scale, 
Lines : integer; CrossHairs : boolean); 


{ Draws a polygon defined in A with "NPoints” points, line style “Line” } 
{ and optional Lines from the axis to the points and cross hairs. } 


procedure Hatch(X_1, Y_1, X_2, Y_2, Delta : real); 
{ Hatch a bar in a histogram } 


procedure DrawHistogram(A :PlotArray; NPoints : integer; 
Hatching : boolean; HatchStyle : integer); 

{ Draw a histogram defined in A with "NPoints" points and optional } 

{ hatching with a given hatch style } 


implementation 


procedure Error{(ProcName : WrkString)}; 
{ Report that an error occurred in the procedure named "ProcName” } 
begin 
DrawString(’/ERROR in ’); 
DrawString(ProcName); 
DrawString(' Press the Button to exit.’); 
repeat until Button; 
Halt; 
end; { Error } 


procedure SetForegroundColor{(Color : integer)}; 
{ Set the foreground drawing color } 
begin 
case Color of 
0 : begin 
ForeColor(BlackColor); 
ForeColorGlb := 0; 


end; 
1 : begin 
ForeColor(WhiteColor); 
ForeColorGlod := 1; 
end; 
end; 


end; { SetForegroundColor } 


procedure SetBackgroundColor{(Color : integer)}; 
{ Set the background color } 
begin 
case Color of 
0 : BackColor(BlackColor); 
1 : BackColor(WhiteColor); 
end; 
end; { SetBackgroundColor } 


procedure DP{(X, Y : integer)}; 
{ Plot a pixel at position (X, Y) } 
begin 

MoveTo(X, Y); 

LineTo(X, Y); 
end; { DP } 


function PD{(X, Y : integer) : boolean}; 
{ Return true if the color of the pixel at (X, Y) matches ForeColorGlb } 
var 
BlackPixel : boolean; 
PixelColor : integer; 
begin 
BlackPixel := GetPixel(X, Y); 
if BlackPixel then 
PixelColor := 0 
else 
PixelColor := 1; 
PD := PixelColor = ForeColorGlb; 
end; { PD } 


procedure DrawStraight{(X1, X2, Y : integer)}; 
{ Draw a horizontal line from X1,Y to X2,Y } 


begin 
MoveTo(X1, Y); 
LineTo(X2, Y); 
end; { DrawStraight } 


procedure InvertWindow{(WinNum : integer)}: 
{ Invert the window referenced by WinNum } 
begin 
if WinNum in [1..MaxWindowGlb] then 
InvertRect(Window[WinNum].W.Port.PortRect) 
else 
Error(’Invert Window’); 
end; { InvertWindow } 


function RealToStr{(R : real) : WrkString}; 
{ Returns the string representation of the real R. } 
var 

Int, Frac : Longint; 

91, S2 : Str255; 

Negitive : boolean; 


begin 
siec= "; 
S2e= "" 
if R < 0.0 then 
Negitive := TRUE 
else 
Negitive := FALSE; 
R := ABS(R); 


Int := Trunc(R); 
Frac := Round(100.0 * (R - Int)); 
NumToString(Int, $1); 
NumToString(Frac, S2); 
if Length(S2) = 1 then 
S2 := S2 + '0'; 
S2 := S1 + '.' + S2; 
if Negitive then 
RealTostr := ’-’ + S2 
else 
RealToStr :=’'’ + S2; 
end; { RealToStr } 


function IntToStr{(l : integer) : WrkString}; 
{ Returns the string representation of the integer I. 
var 
Form : DecForm; 
str : DecStr; 
begin 


Form.Style := FixedDecimal; 
Form.Digits := 0; 
Num2Str(Form, |, Str); 
IntToStr := Str; 

end; { IntToStr } 


procedure SetLineStyle{(Style : integer)}; 
var 
LineStyle : Pattern; 
begin 
case Siyle of 
0 : LineStyle := Black; 
1: LineStyle := White; 
2: LineStyle := Gray; 
3 : LineStyle := LtGray; 
4 : LineStyle := DkGray; 


otherwise 
Style := 0; 
LineStyle := Black; 
end; 
LineStyleGlb := Style; 
PenPat(LineStyle); 


end: { SetLineStyle } 


function GeiLineStyle{ : integer}; 
begin 

GetLinestyle := LineStyleGlb; 
end; { GetLineStyle } 


procedure DefineWorld{(WorldNum : integer; X_1, Y_1, X_2, Y_2 : real)}; 
begin 
if ((X_1 <> X_2) and (Y_1 <> Y_2)) and 
(WorldNum in [1..MaxWorldsGlb]) then 
with World[WorldNum] do 


begin 
X1 := X_1; 
We = Yo 
X2 := X_2;} 
V2 t=) 2; 


if WorldNum > MaxWorldGlb then 
MaxWorldGlb := WorldNum; 

end 

else if WorldNum in [1..MaxWorldsGlb] then 
Error(‘DefineWorld #1’) 

else 
Error(’DefineWorld #2’); 

end; { DefineWorld } 
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procedure SelectWorld{(WorldNum : integer)}; 
begin 
if (WorldNum in [1..MaxWorldGlb]) then 
with World[WorldNum] do 


begin 
WorldNdxGlib := WorldNum; 
X1WidGib := X1; 
Y1WIdGlb := Y1; 
X2WIdGlIb := X2; 
Y2WIdGIb := Y2; 

end 

else 

Error(‘SelectWorld’); 


end; { SelectWorld } 


procedure DefineWindow{(WinNum, XLo, YLo, XHi, YHi, WindowType : integer; Visible : 

boolean)}; 

var 
BoundsRect : Rect; 
Title : Str255; 
RefCon — : Longint; 
Visible : boolean; 
GoAwayFlag : boolean; 


begin 
if WinNum in [1..MaxWindowsGlb] then 
begin 
if WinNum > MaxWindowGlb then 
MaxWindowGlb := WinNum; 
with BoundsRect do 


begin 
Left := XLo; 
Top := YLo; 
Right := XHi; 
Bottom := YHi; 
end; 


Window[WinNum].X1 := XLo; 
Window[WinNum].Y1 := YLo; 
Window[WinNum].X2 := XHi; 
Window[WinNum].Y2 := YHi; 


Title := "3 
GoAwayFlag := TRUE; 
Visible := FALSE; 


RefCon := WinNum; 
if WindowType = documentProc then 
WindowType := 16 * WindowType + 8; { Add Zoom box } 


ie 


Window[WinNum].P := NewWindow(@Window[WinNum].W, BoundsRect, 
Title, Visible, WindowType, 
POINTER(-1), GoAwayFlag, RefCon); 
end 
else 
Error(‘DefineWindow)); 
end; { DefineWindow } 


procedure SetVisibility{((WinNum : integer; Visible : boolean)}; 
begin 

ShowHide(Window[WinNum].P, Visible); 
end; { SetVisibility } 


procedure SelectWind{(WinNum : integer; Visible : boolean)}; 
begin 
if (WinNum in [1..MaxWindowGlb]) then 
begin 
SelectWindow(Window[WinNum].P); 
if Visible then 
ShowWindow(Window[WinNun].P); 
SetPort(@Window[WinNum].W.Port); 
with Window[WinNum] do 
begin 
WindowNdxGib := WinNum; 
X1RefGid := W.Port.PortRect.Left; 
YiRefGid := W.Port.PortReéct. Top; 
X2RefGlib := W.Port.PortRect.Right; 
Y2RefGib := W.Port.PortRect.Bottom; 
BxGlb := (X2RefGlb - X1RefGlb - 16) / (X2WidGlIb - X1WIdGIb); 
ByGlib := (Y2RefGib - Y1RefGlb - 16) / (Y2WIdGIb - Y1WIdGIb); 
AxGlb := X1RefGlb - X1WIidGIb * BxGlb; 
AyGlb := Y1RefGlb - Y1WidGlb * ByGlb; 
AxisGlb := FALSE; 
end; 
end 
else 
Error(‘SelectWind’); 
end; { SeleciWind } 


procedure DefineHeader{(WinNum : integer; Hdr : WrkString)}; 
begin 
Window[WinNum].Header := Hdr; 
SetWTitle(Window[WinNum].P, Hdr); 
end; { DefineHeader } 


procedure RemoveHeader{(WinNum : integer)}; 
begin 
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DefineHKeader(WinNum, "); 
end; { RemoveHeader } 


procedure ReDefineWindow{(WinNum, X1, Y1, X2, Y2, WindowType : integer)}; 
begin 
if (WinNum in [1..MaxWindowsGlb]) then 
begin 
CloseWindow(Window[WinNum].P); 
Window[WinNum].P := NIL; 
DefineWindow(WinNum, X1, Y1, X2, Y2, WindowType); 
SelectWind(WinNum, FALSE); 
DefineHeader(WinNum, Window(WinNum].Header); 
end 
else 
Error('ReDefineWindow'); 
end; { ReDefineWindow } 


function WindowX{(X : real) : integer}; 
var 
Temp : real; 
begin 
Temp := AxGlb + BxGlb * X; 
if Temp > Maxint then 
WindowX := MaxInt 
else if Temp < -32767 then 
WindowX := -32767 
else 
WindowX := trunc(Temp); 
end; { WindowxX } 


function WindowY{(Y : real) : integer}; 
var 
Temp : real; 
begin 
Temp := AyGlb + ByGlb * Y; 
if Temp > Maxlint then 
WindowY := Maxint 
else if Temp < -32767 then 
WindowY := -32767 
else 
WindowY := trunc(Temp); 
end; { WindowY } 


procedure ResetWorlds; 
var 

| : integer; 
begin 


for | := 1 to MaxWorldsGlb do 
DefineWorld(|, XMinGlb, YMinGlb, XMaxGlb, YMaxGlb); 
SelectWorld(1); 
end; { ResetWorlds } 


procedure InitGraphic; 


var 
Index : Integer; 
BoundsFect : Rect; 
Title > Str255: 


Visible : S00lean; 
RefCon —_: Longint; 
GoAwayFlag : boolean; 
WindowType : integer; 
begin 
XMinGlb := screenBits.bounds.Left; 
YMinGlb := screenBits.bounds. Top; 
XMaxGlb := screenBits.bounds.Right; 
YMaxGlb := screenBits.bounds.Bottom; 
for Index := 1 to MaxWindowsGlb do 
begin 
Window([Index].P := NIL; 
Window([index].H] i= NiL; 
end; 
ReseitWorids:; 
MaxWorldGlb := 0; 
MaxWindowGlb := 0; 
WindowNdxGlb := 0; 
WorldNdxGlb := 0; 
AspectFactor := 0.44; 
AspectGlb := ABS(AspectFactor) * AspectFactor; 
AxisGlb := false; 
HatchGlb := false; 
ClippingGlb := true; 
SetLineStyle(0); { Solid Black lines } 
end; { inttGrapnic } 


function Cilip{(var X1, Y1, X2, ¥Y2 : integer) : boolean}; 
var 
Ix1, iy1, Ix2, ly2, Dummy, X1Loc, X2Loc : integer; 
ClipLoc : boolean; 
Temp : real; 


function Inside(X, Xx1, Xx2 : integer) : integer; 
begin 
Inside := 0; 


if X < Xx1 then 


Inside := -1 
else if X > Xx2 then 
Inside := 1; 


end; { Inside } 


begin { Clip } 
Clip := true; 
ClipLoc := true; 
if ClippingGlb then 
begin 
X1Loc := X1RefGlb; 
X2Loc := X2RefGlb; 
Ix1 := Inside(X1, X1Loc, X2Loc); 
lyi := Inside(Y1, Y1RefGlb, Y2RefGlb); 
Ix2 := Inside(X2, X1Loc, X2Loc); 
ly2 := Inside(Y2, Y1RefGlb, Y2RefGlb); 
if (Ix1 or Ix2 or ly1 or ly2) <> 0 then 
begin 
if X1 <> X2 then 
begin 
if x1 <>0O then 
begin 
if Ix1 <0 then 
Dummy := X1Loc 
else 
Dummy := X2Loc; 
if Y2 <> Y1 then 
begin 
Temp := (Y2 - Y1) / (X2 - X1) * (Dummy - X1); 
if Temp > Maxint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp ‘= -32767; 
Y1 := Y1 + trunc(Temp); 
end; 
X1 := Dummy; 
end; 
if (Ix2 <> 0) and (X1 <> X2) then 
begin 
if Ix2 < 0 then 
Dummy := X1Loc 
else 
Dummy := X2Loc; 
if Y2 <> Y1 then 
begin 
Temp := (Y2 - Y1) / (X2 - X1) * (Dummy - X1); 
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if Temp > Maxint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp := -32767; 
Y2 := Y1 + trunc(Temp); 
end; 
X2 := Dummy; 
end; 
ly1 := Inside(Y1, Y1RefGlb, Y2RefGlb); 
ly2 := Inside(¥2, Y1RefGlb, Y2RefGlb); 
end; 
if Yi <> Y2 then 
begin 
if ly1 <> 0 then 
begin 
if ly1 <0 then 
Dummy := Y1RefGlb 
else 
Dummy := Y2RefGlb; 
if X1 <> X2 then 
begin 
Temp := (X2 - X1) / (Y2 - Y1) * (Dummy - Y1); 
if Temp > Maxint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp ‘= -32767; 
X1 := X1 + trunc(Temp); 


end; 
Y1 := Dummy; 
end; 
if ly2 <> 0 then 
begin 
if ly2 < O then 
Dummy := Y1RefGlib 
else 


Dummy := Y2RefGib; 
if Xi <> X2 then 
begin 
Temp := (X2 - X1) / (Y2 - Y1) * (Dummy - Y1); 
if Temp > Maxtint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp := -32767; 
X2 := X1 + trunc(Temp); 
end; 
Y2 := Dummy; 
end; 
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end; 
ly1 := Inside(Y1, Y1RefGlb, Y2RefGlb); 
ly2 := Inside(Y2, Y1RefGlb, Y2RefGlb); 
if (ly1 <> 0) or (ly2 <> 0) then 


ClipLoc := false; 
if ClipLoc then 
begin 


Ix1 := Inside(X1, X1Loc, X2Loc); 
Ix2 := Inside(X2, X1Loc, X2Loc); 
if (Ix2 <> 0) or (Ix1 <> 0) then 


ClipLoc := false; 
end; 
Clip := ClipLoc; 
end; 
end; 
end; { Clip } 
procedure DrawPoint{(Xr, Yr : real)}; 
var 
X, Y : integer; 

begin 


X := WindowX(Xr); 
Y := WindowyY(Yr); 
DP(X, Y); 

end; { DrawPoint } 


function PointDrawn{(Xr, Yr : real) : boolean}; 
begin 

PointDrawn := PD(WindowX(Xr), WindowY(Yr)); 
end; { PointDrawn } 


procedure DrawLine{(X1, Y1, X2, Y2 : real)}; 
begin 
MoveTo(WindowX(X1), WindowY(Y1)); 
LineTo(WindowX(X2), WindowY(Y2)); 
end; { DrawLine } 


procedure DrawLineClipped{(X1, Y1, X2, Y2 : integer)}; 
begin 
if Clip(X1, Y1, X2, Y2) then 
begin 
MoveTo(X1, Y1); 
LineTo(X2, Y2); 
{end 
else 
begin 
Moveto(X1 - 3,Y1 - 3); 
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TextSize(9); 


DrawString(Ds); 
LabelSet := False;} 
end; 


end; { DrawLineClipped } 


procedure DrawCrossDiag{(X, Y, Scale : integer)}; 
begin 
DrawLineClipped(X - Scale, Y + Scale, X + Scale + 1, Y - Scale - 1); 
DrawLineClipped(xX - Scale, Y - Scale, X + Scale + 1, Y + Scale + 1); 
end; { DrawCrossDiag } 


procedure DrawWye{(X, Y, Scale : integer)}; 

begin 
DrawLineClipped(X - Scale, Y - Scale, X, Y); 
DrawLineClipped(X + Scale, Y - Scale, X, Y); 
DrawLineClipped(X, Y, X, Y + Scale); 

end; { DrawWye } 


procedure DrawDiamond{(X, Y, Scale : integer)}; 

begin 
DrawLineClipped(X - Scale, Y, X, Y - Scale - 1); 
DrawLineClipped(X, Y - Scale + 1, X + Scale, Y + 1); 
DrawLineCiipped(X + Scale, Y + 1, X, Y + Scale); 
DrawLineClipped(X, Y + Scale, X - Scale, Y); 

end; { DrawDiamond } 


procedure DrawCircleDirect{(Xr, Yr, R : integer)}; 


type 
Circ = array[1..14] of integer; 


var 
Xk1, Xk2, Yk1, Yk2, Xp1, Yp1, Xp2, Yp2 : integer; 
Xfact, Yfact : real; 


| : integer; 
X : Circ; 
procedure InitX; 
begin 
Al) =e 
X(2] := ar 
Alsi c="23c, 
X(4)=e355- 
X[5] := 465; 
X[G]i=s 56s. 
X[7] := 663; 
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X[8] := 749; 
X[9] := 823; 
X{10] := 885; 
X(11] := 935; 
X[12] := 971; 
X[13] := 993; 
X[14] := 1000; 
end; { InitX } 


begin { DrawCircleDirect } 
InitX; 
Xfact := abs(R); 
Yfact := Xfact * AspectGlb; 
if Xfact > 0.0 then 


begin 
Xk1 := trunc(X[1] * Xfact + 0.5); 
Yk1 := trunc(X[{14] * Yfact + 0.5); 


for | := 2 to 14 do 


begin 
Xk2 := trunc(X{l] * Xfact + 0.5); 
Yk2 := trunc(X{14 - | + 1] * Yfact + 0.5); 
Xp1 := Xr - Xk1; 
Yp1 := Yr + Yk1; 
Xp2 := Xr - Xk2; 
Yp2 := Yr + Yk2; 
DrawLine(Xp1, Yp1, Xp2, Yp2); 
Xp1 i= Xr + Xk1; 
Xp2 := Xr + Xk2; 
DrawLine(Xp1, Yp1, Xp2, Yp2); 
od c= Yr - YK1; 
Wed := Yr - Yk2; 
DrawLine(Xp1, Yp1 + 1, Xp2, Yp2 + 1); 
Xp1 := Xr - Xk1; 
Xp2 := Xr - Xk2; 
DrawLine(Xp1, Yp1 + 1, Xp2, Yp2 + 1); 
Xk1 := Xk2; 
Yk1 := Yk; 
end; 
end 
else 
mrixXr, Yr); 


end; { DrawCircleDirect } 


procedure DrawCircle{(x_R, Y_R, Xradius : real)}; 
begin 

DrawCircleDirect(WindowX(X_R), WindowY(Y_R), trunc(Xradius)); 
end; { DrawCircle } 
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procedure DrawCross{(X1, Y1, Scale : integer)}; 
begin 
DrawLineClipped(X1 - Scale, Y1, X1 + Scale + 2, Y1); 
DrawLineClipped(X1, Y1 - Scale, X1, Y1°+ Scale + 1): 
end; { DrawCross } 


procedure DrawStar{(X, Y, Scale : integer)}; 

begin 
DrawLineClipped(X - Scale, Y + Scale, X + Scale + 1, Y - Scale - 1); 
DrawLineClioped(X - Scale, Y - Scale, X + Scale + 1, Y + Scale + 1); 
DrawLineClipped(X - Scale - 2, Y, X + Scale + 4, Y); 

end; { DrawStar } 


procedure DrawSquareC{(X1, Y1, X2, Y2 : integer; Fill : boolean)}; 
var 
| : integer; 


procedure DSC(X1, X2, Y : integer); 
begin 

DrawStraight(X1, X2, Y); 
end; { DSC } 


begin { DrawSquareC } 
if not Fill then 
begin 
DrawLineClipped(X1, Y1, X2, Y1 


( i: 

DrawLineClipped(X2, Y1, X2, Y2); 

DrawLineClipped(X1, Y2, X2, Y2); 
( ) 


DrawLineClipped(X1, Y2, X1, Y1 


’ 


end 
else 
for | := Y2 to Y1 do 
DSC(X1, X2, |); 


end; { DrawSquareC } 


procedure DrawSquare{(X1, Y1, X2, Y2 : real; Fill : boolean)}; 
var 

|, X1Loc, YiLoc, X2Loc, Y2Loc : integer; 

DirectModeLoc : boolean; 


procedure DS(X1, X2, Y : integer); 
begin 
if LineStyleGlb = O then 
DrawStraight(X1, X2, Y) 
else 
DrawLine(X1, Y, X2, Y); 
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end; { DS } 


procedure DSC(X1, X2, Y : integer); 
begin 

DS(X1, X2, Y); 
end; { DSC } 


procedure DrawSqr(X1, Y1, X2, Y2 : integer; Fill : boolean); 
var 
| : integer; 
begin 
if not Fill then 
begin 
DS(X1, X2, Y1); 
DrawLine(X2, Y1, X2, Y2); 
DS(X1, X2, Y2); 
DrawLine(X1, Y2, X1, Y1); 


end 
else 
for | := Y1 to Y2 do 
DS(X1, X2, |); 


end; { DrawSar } 


begin { DrawSquare } 
XiLoc := WindowX(X1); 
Y1Loc := WindowY(Y1); 
X2Loc := WindowX(X2); 
Y2Loc := WindowY(Y2) 
if not Fill then 
begin 
DSC(X1Loc, X2Loc, Y1Loc); 
DrawLineClipped(X2Loc, Y1Loc, X2Loc, Y2Loc); 
DSC(X1Loc, X2Loc, Y2Loc); 
DrawLineClipped(X1Loc, Y2Loc, X1iLoc, Y1Loc); 
end 
else 
for | := Y1Loc to Y2Loc do 
DSC(X1Loc, X2Loc, |); 
end; { DrawSquare } 


’ 


procedure DrawAscii{(X, Y : integer; Size, CharByte : byte)}; 
begin 
MoveTo(X, Y); 
TextSize(Size * 12); 
DrawChar(Chr(CharByte)); 
end; { DrawAscii } 
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procedure DrawText{(X, Y, Scale : integer; Txt : WrkString)}; 
var 
Index : integer; 
EscStr : boolean; 
StringLen : integer; 
AsciiValue : integer; 
SymbolScale : integer; 
SsymbolCode : integer; 
begin 
Index := 1; 
EscSir i= FALSE; 
StringLen := Length(Txt); 
while (Index <= StringLen) and (not EscStr) do 
begin 
if Txt[Index] = #27 then 
EscStr := TRUE; 
Index := Index + 1; 
end; 
if not EscStr then 
begin 
MoveTo(X, Y); 
TextSize(Scale * 12); 


DrawString(T xt); 
end 
else 
begin 
Index := 1; 
while Index <= StringLen do 
begin 


AsciiValue := Ord(Txt[Index]); 
if AsciiValue = 27 then 
begin 
SymbolScale := Scale; 
Index := Index + 1; 
if Index <= StringLen then 
begin 
SymbolCode := Ord(Txt[Index]) - 48; 
if (Index + 2 <= StringLen) and (Ord(Txt[Index + 1]) = 64) then 
begin 
SymbolCode := Ord(Txt[Index]) - 48; 
Index := Index + 2; 
end; 
case SymbolCode of 
1 : DrawCross(X + SymbolScale, Y + Scale, SymbolScale); 
2 : DrawCrossDiag(X + SymbolScale, Y + Scale, SymbolScale); 
3,4 : DrawSquareC(X, Y + (SymbolScale shi 1) - 1, 
X + (SymbolScale shl 1), Y - 1, (SymbolCode = 4)); 
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5 : begin 
DrawDiamond(X + trunc(1.5 * SymbolScale), 
Y + SymbolScale - 1, SymbolScale + 1); 
X := X + SymbolScale; 
end; 
6 :DrawWye(X + SymbolScale, Y + SymbolScale - 1, SymbolScale); 
7 : begin 
DrawStar(X + SymbolScale shi 1, Y + SymbolScale - 1, SymbolScale); 
X := X + SymbolScale shl 1; 
end; 
8 : DrawCircleDirect(X + SymbolScale, Y + (SymbolScale shr 1), 
SymbolScale + 1); 
end; 
X := X + 3 * SymbolScale; 
SymbolScale := Scale; 
end; 
end 
else 
DrawAscii(X, Y, Scale, AsciiValue); 
Index := Index + 1; 
end; 
end; 
end; { DrawText } 


procedure DrawTextW{(X, Y : real; Scale : integer; Txt : WrkString)}: 
begin 

DrawText(WindowX(X), WindowY(Y), Scale, Txt); 
end; { DrawTextW } 


procedure TextStyle{(Face : Style)}; 
{ Face = (bold, italic, underline, outline, shadow, condense, extend) } 
begin 
TextFace(Face); 
end; { TextStyle } 


procedure HardCopy{(TopWin : boolean)}; 
begin 
PrDrvrOpen; 
if TopWin then 
{ Print the top folder. } 
PrctiCall(iPrEvtCtl, LPrEvtTop, 0, LScreenBits) 
else 
{ Print the whole screen. } 
PrctiCall(iPrEvtCtl, LPrEvtAll, 0, LScreenBits); 
PrDrvrClose; 
end; { HardCopy } 
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procedure OpenPic{(WinNum : integer; ShowPic : boolean)}; 
begin 
if Window{WinNum].H <> NIL then 
begin 
KillPicture(Window[WinNum].H); 
Window[WinNum].H := NIL; 
end; 
RectRgn(Window[WinNum].W.Port.clipRgn, ScreenBits.bounds); 
Window[WinNum].H := OpenPicture(Window[WinNum].W.Port.PortRect); 
if ShowPic then 
ShowPen 
end; { OpenPic } 


procedure DrawPic{(WinNum : integer)}; 
var 
PictRect : Rect; 
begin 
PictRect := Window[WinNum].W.Port.PortRect; 
with PictRect do 
begin 
if ((Bottom - Top) < 200) OR ((Right - Left) < 200) then 
begin 
Right := Right - 16; { so we don't overwrite the grow region } 
Bottom := Bottom - 16;{ on a small window. } 
end; 
end; 
DrawPicture(Window[WinNum].H, PictRect); 
end; { DrawPic } 


procedure ErasePic{(WinNum : integer)}; 
begin 

if Window[WinNum].H <> NIL then 

begin 

KillPicture(Window[WinNum].H); 
Window[WinNum].H := NIL; 

end; 

end; { ErasePic } 


procedure ClearWindow{(WinNum : integer)}; 
begin 

EraseRect(Window[WinNum].W.Port.PortRect); 
end; { ClearWindow } 


function WhereX{ : integer}; 
var 

Pt : Point; 
begin 
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GetPen(Pt); 
WhereX := Pt.H; 
end; { WherexX } 


function WhereY{ : integer}; 
var 
Pt : Point; 
begin 
GetPen(Pt); 
WhereY := Pt.V; 
end; { WhereY } 


procedure SetWindow{(X1, Y1, X2, Y2 : integer)}; 


begin 
X1RefGlb := X1; 
YiRefGlb := Y1; 
X2RefGlb := X2; 
Y2RefGlb := Y2; 
BxGlb := (X2 - X1) / (X2WIdGIb - X1WlidGlb); 
ByGlb := (Y2 - Y1) / (Y2WIdGIb - Y1WidGlb); 


AxGlb := X1 - X1WidGlb * BxGlb; 
AyGlb := Y1 - Y1WidGlb * ByGlb; 
AxisGlb := FALSE; 

end; { SetWindow } 


procedure FindWorld{(| : integer; A : PlotArray; NPoints : integer)}; 


var 
J : integer; 
Xmax, Ymax, Xmin, Ymin, Xmid, Ymid, Xdiff, Ydiff : real; 


begin 
NPoints := abs(NPoints); 
if NPoints > 2 then 
if | in [1..MaxWorldsGlb] then 


begin 
Xmax := A[t, 1]; 
Ymax := A[1, 2]; 


Xmin := Xmax; 
Ymin := Ymax; 
for J := 2 to NPoints do 
begin 
if A[J, 1] > Xmax then 
Xmax := Af, 1] 
else 
if A[J, 1] < Xmin then 
Xmin := AlJ, 1]; 
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if A[J, 2] > Ymax then 


Ymax := Alfu, 2] 
else 
if A{J, 2} < Ymin then 
Ymin := Afu, 2]; 
end; 


Xmin := Round(Xmin); 

Ymin := Round(Ymin) - 0.5; 
Xmax := Round(Xmax); 
Ymax := Round(Ymax) + 0.5; 

DefineWor'd(i, Xmin, Ymin, Xmax, Ymax); 

SelectWorld(l); 

end 
else 

Error(’FindWorld #1’) 

else 
Error(’FindWorld # 2’); 
end; { FindWorld } 


procedure FindWorld1{(I : integer; A : PlotArray; NPoints : integer)}; 


var 
J : integer; 
Xmax, Ymax, Xmin, Ymin, Xmid, Ymid, Xdiff, Ydiff : real; 


begin 
Xmin := XMn; 
Ymin := YMn{Round(YMn)} ; 
Xmax := XMx{Round(XMx)}; 
Ymax := YMx{Round(YMx)} ; 
DefineWorld(|, Xmin, Ymin, Xmax, Ymax); 
SelectWorld(l); 
end; { FindWorld } 


procedure DrawAxis{(Footer1, Footer2 : WrkString; Arrows : boolean)}; 
var 

LineStyleLoc, XkO, YkO, Xk1, Yk1, Xk2, Yk2, 

MaxExponentX, MaxExponentyY, TickPoint : integer; 

TickSmali, TickLarge, Max, Min, Tick, Offset, Diff : real; 


function Log(X : real) : real; 
{ Base 10 logarithm of X. } 


begin 
Log := Ln(X) / Ln(10.0); 
end; { Log } 


function ALog(X : real) : real; 
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{ Ten raised to the X power. } 


begin 
ALog := Exp(X * Ln(10.0)); 
end; { Alog } 


function Frac(R : real) : real; 
{ Return the fractional part of the real number R. } 


begin 
Frac := R - Int(R); 
end; { Frac } 


procedure Ticks(NTicks : integer; Max, Min : real; 
var TickSmall, TickLarge : real); 

{ NTicks : The approximate number of tick marks in the interval. } 
{ Max, Min : World coordinates of the axis extremes. } 
{ TickSmall, TickLarge :Tick mark intervals, in world coordinates .} 
var 

TickLog : array[1..4] of real; 

|, J, ChA : integer; 

Delta, XTicks, LogTicks, Mant, MinDiff, Diff : real; 


begin 
TickLog[1] := 0.0; 
TickLog[2] := Log(2.0); 
TickLog[3] := Log(5.0); 


TickLog[4] := 1.0; 
XTicks := NTicks; 
Delta := Max - Min; 
XTicks := Delta / XTicks; 
LogTicks := Log(XTicks); 
ChA := trunc(LogTicks); 
if LogTicks < 0.0 then 
ChA := ChA - 1; 
MinDiff := 1.0; 
Mant := LogTicks - ChA; { Fractional part of logarithm } 
for | := 1 to 4 do 
begin 
Diff := Abs(Mant - TickLog[l]); 
if (Diff < MinDiff) then 
begin 
MinDiff := Diff; 
J := I; 
end; 
end; 
LogTicks := ChA + TickLog[J]; { Logarithm of tick mark } 
TickLarge := ALOg(LogTicks); { The tick mark } 


{ Find the small tick marks, that are two tick scales smaller } 
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if J > 2 then 


begin 
J She 2; 
LogTicks := ChA + TickLog|[J}; 
end 
else 
begin 
J:=J+1; 
LogTicks := ChA + TickLog[J] - 1; 
end; 
TickSmail := ALog(LogTicks); 
end; { Ticks } 


function StringNumber(X1 : real; MaxExponent : integer) : WrkString; 
begin 

StringNumber := RealToStr(X1 * Exp(-MaxExponent * Ln(10.0))); 
end; { StringNumber } 


function GetExponent(X1 : real) : integer; 
begin 
GetExponent := 0; 
if X1 <> 0.0 then 
if ABS(X1) >= 1.0 then 
GetExponent := trunc(Ln(ABS(X1)) / Ln(10.0)) 
else 
GetExponent := -trunc(ABS(Ln(ABS(X1))) / Ln(10.0) + 1.0); 
end; { GetExponent } 


procedure DrawNum(X1, Y1, MaxExponent : integer; Number : real); 
var 
StrNumber : WrkString; 
begin 
TextSize(9); 
StrNumber := StringNumber(Number, MaxExponent); 
Cee Ol = 6 
MoveTo(X1, Y1); 
DrawString(StrNumber); 
TextSize(12); 
end; { DrawNum } 


procedure DrawExponent(X1, Yi, MaxExponent : integer); 
var 

NumStr : WrkString; 
begin 

MoveTo(X1, Y1); 

TextSize(9); 

DrawChar('x’'); 
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DrawChar(’ '); 
DrawChar('1’); 
DrawChar(’0’); 
X1 := WherexX + 1; 
Y1 := WhereY - 3; 
MoveTo(X1, Y1); 
NumStr := IntToStr(MaxExponent); 
TextSize(7); 
DrawString(NumStr); 
TextSize(12); 
end; { DrawExponent } 


begin { DrawAxis } 
LineStyleLoc := LinestyleGlb; 
SetLineStyle(0); { Black } 


XkO := X1RefGib + X1 Offset; 
YkO := Y2RefGlb - Y2Offset; 
Xk1 := XkO; 
Yk1 := Y1RefGib + Y1Offset; 
Xk2 := X2RefGib - X2Offset; 
Yk2 := YkO; 


MoveTo(XKO, YKO); { Draw the Y axis with optional Arrows } 
LineTo(XK1, YK1); 
if Arrows then 
begin 
MoveTo(Xk0, Yk1); 
LineTo(XkO - 4, Yk1 + 4); 
MoveTo(Xk0, Yk1); 
LineTo(XkO + 4, Yk1 + 4); 
DP(XkO, Yk1 - 1); 
end; 


MoveTo(Xk0O, YkOQ); { Draw the X axis with optional Arrows } 
LineTo(Xk2 + 1, Yk2); 
if Arrows then 
begin 
MoveTo(Xk2, Yk2); 
LineTo(Xk2 - 4, Yk2 - 4); 
MoveTo(Xk2, Yk2); 
LineTo(Xk2 - 4, Yk2 + 4); 
end; 


if Footer! <> “ then { Draw the 1st footer below the X axis } 


begin 
MoveTo(Xk0, YkO + 45); 
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TextSize(9); 
DrawString(Footer1); 
end; 


if Footer2 <> ” then { Draw the 2nd footer below the X axis } 
begin 

MoveTo(Xk0, YkO + 65); 

TextSize(9); 

DrawString(Footer2); 
end; 


if (ABS(YkO - Yk1) >= 35) and (ABS(Xk2 - Xk1) >= 150) then 
begin 
if ABS(Y2WIdGIb) > ABS(Y1WIdGlb) then 
MaxExponentY := GetExponent(Y2WIdGlb) 
else 
MaxExponentyY := GetExponent(Y1WIdGlb); 


if MaxExponentY <> 0 then { Draw the power of ten on top of Y axis } 
DrawExponent(Xk1 - 30, Yk1 + 2, MaxExponentyY); 


TickPoint := YkO; 
if YIWIidGlb > Y2WidGib then 
begin 
Max := YIWIdGib; 
Min := Y2WidGie; 
end 
else 
begin 
Max := Y2WIdGlb; 
Min := Y1WidGlb; 
end; 


{ Using the Max and Min values, this procedure call calculates } 
{ large and small Tick Marks in world coordinates. } 
Ticks(5, Max, Min, TickSmall, TickLarge); 


Offset := Min / TickLarge; 
Offset := Offset - Frac(Offset); 
Tick := Offset * TickLarge; 
if Tick < Min then 

Tick := Tick + TickLarge; 


{ Tick is the world coordinate at which the tick mark is to be drawn } 
Diff := Max - Min; 

{ Plot large tick marks and Numeric labels } 

while Tick <= Max do 
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begin 
TickPoint := YkO - Trunc((YkO - Yk1) * (Tick - Min) / Diff); 
MoveTo(XkO, TickPoint); 
LineTo(XkO - 4, TickPoint); 
DrawNum(X1Offset - 30, TickPoint + 7, MaxExponentyY, Tick); 
Tick := Tick + TickLarge; 
end; 


{ The same repeated for the small tick marks, } 
{ only without axis numbering. } 
Offset := Min / TickSmall; 
Offset := Offset - Frac(Offset); 
Tick := Offset * TickSmall; 
if Tick < Min then 
Tick := Tick + TickSmall; 
while (Tick + 0.01) < Max do 
begin 
TickPoint := YkO - Trunc((YkO - Yk1) * (Tick - Min) / Diff); 
MoveTo(XkO, TickPoint); 
LineTo(XkO - 2, TickPoint); 
Tick := Tick + TickSmall; 
end; 


if ABS(X2WIdGIb) > ABS(X1WldGlb) then 
MaxExponentX := GetExponent(X2WIdGIb) 
else 
MaxExponentX := GetExponent(X1WIdGIb); 
if MaxExponentX <> 0 then { Draw power of ten label on X axis } 
DrawExponent(Xk2 - 25, YkO + 28, MaxExponentx); 


{ This is the same as for the Y axis, but the window } 
{ and world are appropriate for the X axis. } 
TickPoint := Xk0; 
if X1WIdGlb > X2WidGlb then 
begin 
Max := X1WIdGIb; 
Min := X2WIdGlb; 
end 
else 
begin 
Max := X2WIdGIb; 
Min := X1WidGlb; 
end; 
Ticks(5, Max, Min, TickSmall, TickLarge); 
Offset := Min / TickLarge; 
Offset := Offset - Frac(Offset); 
Tick := Offset * TickLarge; 
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if Tick < Min then 
Tick := Tick + TickLarge; 
Diff := Max - Min; 
while Tick <= Max do 
begin 
TickPoint := XkO + Trunc((Xk2 - XkQ) * (Tick - Min) / Diff); 
MoveTo(TickPoint, YkO); 
LineTo(TickPoint, YkO + 4); 
DrawNum(TickPoint - 14, YkO + 20, MaxExponentx, Tick); 
Tick t= Tick + TickLarge; 
ena; 
Offset := Min / TickSmall; 
Offset := Offset - Frac(Offset); 
Tick := Offset * TickSmall; 
if Tick < Min then 
Tick := Tick + TickSmall; 
while (Tick + 0.01) < Max do 
begin 
TickPoint := XkO + Trunc((Xk2 - Xk0) * (Tick - Min) / Diff); 
MoveTo(TickPoint, YkO); 
LineTo(TickPoint, YkKO + 2); 
Tick := Tick + TickSmall; 
end; 
end; 
SetLineStyle(LineStyleLoc); 
AxisGibd := TRUE; 
end; { DrawAxis } 


procedure ResetAxis; 
begin 

AxisGlb := true; 
end; { ResetAxis } 


procedure DrawPolygon{(A : PlotArray; First, NPoints, Line, Scale, 
Lines : integer; CrossHairs : boolean)}; 
var 
i, X1, X2, Yi, Y2, XOffset, YOffset, 
X1Reiloc, YiRefloc, X2RefLoc, Y2RefLoc, 
DeltaY, XOs1, XOs2, YOsi, YOs2 : integer; 
AutoClip, DirectModeLoc, PlotLine, PlotSymbol, Flipped : boolean; 
XiLoc, YiLoc, X2Loc, Y2Loc : integer; 
Temp : real; 
LineStyleLoc2 : integer; 
DrawPt : Boolean; 


procedure DrawPointClipped(X, Y : integer); 
begin 
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if (X1 > X1RefGlb) and (X2 < X2RefGlb) then 
if (Y1 > Y1RefGlb) and (Y2 < Y2RefGlb) then 
DP(X, Y); 
end; { DrawPointClipped } 


procedure Drawltem(X, Y : integer); 
var 
LineStyleLoc : integer; 
begin 
LineStyleLoc := LineStyleGlb; 
SetLineStyle(0); { Black } 
case Line of 
2  : DrawCrossDiag(X, Y, Scale); 
3, 4 : DrawSquareC(X - Scale, Y + Scale, X + Scale, Y - Scale, (Line = 4)); 
5  : DrawDiamond(X, Y, Scale + 1); 
6 : DrawWye(X, Y, Scale + 1); 
1  :DrawCross(X, Y, Scale); 
8  : DrawCircleDirect(X, Y, Scale + 1); 
9 : begin 
PlotLine := false; 
if AutoClip then 
DrawPointClipped(X, Y) 
else 
DP(X, Y); 
end; 
7 +: DrawPoint(X, Y){, Scale)}; 
end; 
SetLineStyle(LineStyleLoc); 
end; { Drawitem } 


begin { DrawPolygon } 
if AxisGlb then 
Flipped := FALSE 
else 
begin 
Flipped := TRUE; 
Temp := World[WorldNdxGlb].¥1 ; 
World[WorldNdxGlb].Y1 := World[WorldNdxGib].Y2; 
World[WorldNdxGlb].Y2 := Temp; 
SelectWorld(WorldNdxGlb); 
SelectWind(WindowNdxGlb, TRUE); 
end; 
if abs(NPoints - First) >= 2 then 
begin 
AutoClip := (NPoints < 0); 
NPoints := abs(NPoints); 
KOS|2= 1: 
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XOs2 := 1; 


YOs1 := 6; 
YOs2 := 6; 

if AxisGlb then 
begin 


XOs1 := X1Offset; 
XOs2 := X2Offset; 
YOs1 := Y10Offset; 
YOs2 := Y2Offset; 
if ({(X2RetGlo - XOs2 - X1RefGlb + XOs1) > (XOs1 + XOs2)) and 
((Y2RefGib - YOs2 - Y1RefGlb + YOs1) > (YOs1 + YOs2)) then 
begin 
X1RefLoc := X1RefGlb; 
X1 := X1RefGib + XOs1; 
Y1Refloc := Y1RefGlb; 
Y1 := Y1RefGlb + YOs1; 
X2RefLoc := X2RefGib; 
X2 := X2RefGlb - XOs2; 
Y2RefLoc := Y2RefGlb; 
Y2 := Y2RefGlb - YOs2; 
SetWindow(X1, Y1, X2, Y2); 
AxisGlb := TRUE; 
end: 
enc; 
PiotLine := (Line >= 0); 
PlotSymboi := (Line <> 0); 
Line := abs(Line); 
Scale := abs(Scale); 
if Lines < 0 then 
DeltaY := Trunc(1.0 / (abs(Y1WIdGIb) + abs(Y2WIdGIb)) * 
abs(Y1WIdGlb) * abs(Y2RefGlb - Y1RefGlb)) + 1 
else 
DeltaY := 0; 
if (NPoints < 2) then 
Error(’‘DrawPolygon #1’) 
else 
begin 
if CrossHairs then 
begin 
LineStyleLoc2 := LineStyleGlb; 
setLineStyle(3); { Light Gray } 
MoveTo(X1RefGib, Y1RefGlb + Y2RefGlb - WindowY(0.0)); 
LineTo(X2RefGlb, Y1RefGlb + Y2RefGlb - WindowY(0.0)); 
MoveTo(WindowX(0.0), Y1RefGlb); 
LineTo(WindowX(0.0), Y2RefGlb); 
SetLineStyle(LineStyleLoc2); 
end; 
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X1 := WindowX(A|First, 1]); 
Y1 := Y1RefGlb + Y2RefGlb - WindowY(A[First, 2]) ; 
Drawltem(X1, Y1); 
if Abs(Lines) = 1 then 
if AutoClip then 
DrawLineClipped(X1, Y2RefGlb - DeltaY, X1, Y1) 
else 
begin 
MoveTo(X1, Y2RefGlb - DeltaY); 
LineTo(X1, Y1); 
end; 
DrawPt := True; 
for |:= First + 1 to NPoints do 
begin 
X2 := WindowX{All, 1]); 
Y2 := Y2RefGib + Y1RefGlb - WindowY(A[l, 2]); 
Drawltem(X2, Y2); 
if StepA then begin 
if not AMarkStatus then 


begin 
if (A[I, 2] > 0) and DrawPt 


begin 
if ARJustification then 


begin 
MoveTo(X2+15, Y2+5); 
TextSize(9); 
DrawString(Ds); 
end 
else 
begin 
Move To(X2-65, Y2); 
TextSize(9); 
DrawString(Ds); 
end; 
DrawPt := False; 
end; 
end 
else 
begin 
if (| > (NPoints - InitDegree)) and (A[I, 2] > 0) and DrawPt 


then {Clip(X1,Y1,X2,Y2)} 


then 


begin 
if ARJustification then 


begin 
MoveTo(X2+15, Y2+5); 


TextSize(9); 
DrawString(Ds); 
end 
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else 
begin 
MoveTo(X2-65, Y2); 
TextSize(9); 
DrawString(Ds); 
end; 
DrawPt := False; 
end; 
end; 
end:{StepA} 
if not StepA then begin 
if not BMarkStatus then 
begin 
if (A{l, 2] > 0) and DrawPt then {Clip(X1,Y1,X2,Y2)} 
begin 
if BRJustification then 
begin 
MoveTo(X2+15, Y2+5); 
TextSize(9); 
DrawString(Ds); 
end 
else 
begin 
Move To(X2-65, Y2); 
TextSize(9); 
DrawString(Ds); 
end; 
DrawPt := False; 
end; 
end 
else 
begin 
if (1 > (NPoints - InitDegree)) and (A[I, 2] > 0) and DrawPt then 
begin 
if BRJustification then 
begin 
MoveTo(X2+15, Y2+5); 
TextSize(9); 
DrawString(Ds); 
end 
else 
begin 
MoveTo(X2-65, Y2); 
TextSize(9); 
DrawString(Ds); 
end; 
DrawPt := False; 


end; 
end; 
end;{not StepA} 


if Abs(Lines) = 1 then 
if AutoClip then 
DrawLineClipped(X2, Y2RefGib - DeltaY, X2, Y2) 
else 
begin 
MoveTo(X2, Y2RefGib - DeltaY); 
LineTo(X2, Y2); 
end; 
if PlotLine then 
if AutoClip then 
DrawLineClipped(X1, Y1, X2, Y2) 
else 
begin 
MoveTo(X1, Y1); 
LineTo(X2, Y2); 


end; 
X1 := X2; 
Vi = Y2; 
end; 
end; 
if AxisGlb then 
begin 


SetWindow(X1RefLoc, Y1RefLoc, X2RefLoc, Y2RefLoc); 
AxisGlb := false; 
end; 
end 
else 
Error(DrawPolygon # 1°); 
if Flipped then 
begin 
Temp := World{[WorldNdxGlb].¥1; 
World(WorldNdxGlb].Y1 := World{[WorldNdxGlb].¥Y2; 
World[WorldNdxGlb].Y2 := Temp; 
SelectWorld(WorldNdxGlb); 
SelectWind(WindowNdxGib, TRUE); 
end; 
end; { DrawPolygon } 


procedure Hatch{(X_1, Y_1, X_2, Y_2, Delta : real)}; 


Var 


X1, Y1, X2, Y2 : integer; 


procedure HatchDirect(X1, Y1, X2, Y2, Delta : integer); 
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var 
|, Yst, Yen, Count : integer; 
X1RefLoc, X2RefLoc, YiRefLoc, Y2RefLoc : integer; 
ClippingLoc : boolean; . 
X1D, Y1D, X2D, Y2D : integer; a 


begin { HatchDirect } 
if Delta <> O then 
begin 
HatchGid := true; 
ClippingLoc := ClippingGlb; 
CliopingGlb := true; 
X1RefLoc := X1RefGlb; 


X1RefGlb := X1; 
X2RefLoc := X2RefGlb; 
X2RefGlb := X2; 
Y1RefLoc := Y1RefGlb; 
Y1RefGlb := Y1; 
Y2RefLoc := Y2RefGlb; 
Y2RefGlb := Y2; 


Yst := Y1 + Delta; 

Yen := Y1 - X2 + X1 + Delta; 
if Delta < O then 
begin 

Delta := -Delta; 

f:= Yst; 

Yst := Yen; 

Yen := |; 
end; 

Count := (Y2 - Y1 + X2 - X1 + X2 - X1) div Delta; 
for | := 1 to Count-1 do 


begin 
X1D := X1; 
Y1D := Yst; 
X2D := X2; 
Y2D := Yen; 
if Clip(X1D, Y1D, X2D, Y2D) then 
begin 


VioveTo(X1D, Y1D); 
LineTo(X2D, Y2D); 
end; 
Yst := Yst + Delta; 
Yen := Yen + Delta; 
end; 
ClippingGlb := ClippingLoc; 
HatchGlb := false; 
X1RefGlb := X1RefLoc; 
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X2RefGlb := X2RefLoc; 
Y1RefGlib := Y1RefLoc; 
Y2RefGlb := Y2RefLoc; 
end; 
end; { HatchDirect } 


begin { Hatch } 
HatchDirect(trunc(X_1), trunc(Y_1), trunc(X_2), trunc(Y_2), trunc(Delta)) 
end; { Hatch } 


procedure DrawHistogram{(A :PlotArray; NPoints : integer; 
Hatching : boolean; HatchStyle : integer)}; 


var 
X1, X2, Y2, NPixels, Delta, NDiff, YRef, LineStyleLoc, | : integer; 
Fract, S, Y : real; 
DirectModeLoc, Negative : boolean; 
X1Loc, Y1Loc, X2Loc, Y2Loc : integer; 
X1RefLoc, YiRefLoc, X2RefLoc, Y2RefLoc, YAxis : integer; 
Temp : real; 


begin { DrawHistogram } 
if ABS(NPoints) >= 2 then 
begin 
LineStyleLoc := LinestyleGlb; 
SetLineStyle(0); { Black } 
if AxisGlb then 
begin 
X1RefLoc := Window[WindowNdxGlb}].X1; 
Y1RefLoc := Window[WindowNadxGlb].Y1; 
X2RefLoc := Window[WindowNdxGlb].X2; 
Y2RefLoc := Window[WindowNdxGlb].Y 2; 
SetWindow(X1RefGlb + X1Offset, Y1RefGib + Y1Offset, 
X2RefGlb - X2Offset, Y2RefGlb - Y2Offset); 
AxisGlb := TRUE; 
end; 
Negative := NPoints < 0; 
NPoints := ABS(NPoints); 
NPixels := X2RefGlb - X1RefGlb; 
Delta := NPixels div NPoints; 
NDiff := NPixels - Delta “ NPoints; 
Fract := NDiff / NPoints; 
S := -Fract; 
X1 := X1RefGlb; 
Temp := Y2RefGlb + Y1RefGlb - AyGlb; 
if Temp > Maxint then 
Temp := Maxint 
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else 
if Temp < -32767 then 
Temp := -32767; 
YRef := trunc(Temp); 
if Negative then 
DrawStraight(X1, X2RefGlb, YRef); 
YAxis := Y1RefGlb; 
if BYGIb > 0 then 
YAxis := Y2RefGlb; 
for | := 71 to NPoints do 
begin 
X2 := X1 + Delta; 
Y¥Y = All, 2]; 
if not Negative then 
Y := ABS(Y); 
Temp := AyGlb + ByGlb * Y; 
if Temp > Maxint then 
Temp := Maxint 
else 
if Temp < -32767 then 
Temp := -32767; 
Y2 := Y2RefGib + Y1RefGlb - trunc(Temp); 
if not Negative then 
begin 
MoveTo(X1, YAxis); 
LineTo{xt, Y2); 
MoveTo(X1, Y2); 
LineTo(X2, Y2); 
MoveTo(X2, Y2); 
LineTo(X2, YAxis); 
if Hatching then 
if Odd(l) then 
Hatch(X1, Y2, X2, YAxis, HatchStyle) 
else 
Hatch(X1, Y2, X2, YAxis, -HatchStyle); 
end 
else 
begin 
MoveTo(X1, YRef); 
LineTo(X1, Y2); 
MoveTo(X1, Y2); 
LineTo(X2, Y2); 
MoveTo(X2, Y2); 
LineTo(X2, YRef); 
if Hatching then 
if YRef - Y2 < 0 then 
if Odd(l) then 
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Hatch(X1, YRef, X2, Y2, HatchStyle) 
else 
Hatch(X1, YRef, X2, Y2, -HatchStyle) 
else 
if Odd(l) then 
Hatch(X1, Y2, X2,YRef, HatchStyle) 
else 
Hatch(X1, Y2, X2, YRef, -HatchStyle); 
end; 
X1 i= X2; 
end; 
if AxisGlb then 
begin 
SetWindow(X1RefLoc, Y1RefLoc, X2RefLoc, Y2RefLoc); 
AxisGlb := FALSE; 
end; 
SetLineStyle(LineStyleLoc); 
end 
else 
Error(‘DrawHistogram’); 
end; { DrawHistogram } 


begin 
end. { TurboGraph } 


* This is the resource file that defines the menus and icons for MacRootLocus. 
MacRootLocus.Rsrc 
TYPE DLOG 


256 (36) 
CE Parameter Dialog 
70 100 300 412 
Visible NoGoAway 
1 
0 
256 


TYPE DITL 
256 (36) 
11 


Btnltem Enabled 
185 240 210 300 
Cancel 


Btnltem Enabled 
150 240 175 300 
CK 


EditTextltem Enabled 
60 260 75 280 


EditTextitem Enabled 
95 165 110 205 
1 


EditTextltem Enabled 
120 165 135 205 
100 


EditTextltem Enabled 
145 165 160 205 
1E-6 


StatText Disabled 
20 15 40 250 
Characteristic Equation Parameter 


StatText Disabled 
60 40 75 240 
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Degree of ths polynomial 


StatText Disabled 
95 40 110 150 
InitG@uess 


StatText Disabled 
120 40 135 150 
Maxiter 


StatText Disabled 
145 40 160 150 
Tolerance 


TYPE DLOG 


257 (36) 
One Parameter Root Locus Plot Data 
50 86 315 426 
Visible NoGoAway 
; 
0 
257 


TYPE DITL 

257 (36) 
21 
Btnitem Enabled 
35 270 55 330 
Cancel 


Btnitem Enabled 
10 270 30 330 
PLOT 


Radioltem Enabled 
95 40 110 265 
Linear Point Interval 


Radioltem Enabled 

115 40 130 265 
Logarithmic Point Interval 
Radioltem Enabled 

145 25 160 150 
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Auto Scale Axis 


Radioltem Enabled 
145 165 160 320 
Manual Scale Axis 


EditTextitem Enabled 
65 110 80 160 
0.1 


EditTextitem Enabled 
65 265 80 315 
10000 


EditTextltem Enabled 
170 85 185 125 
-10 


EditTextltem Enabled 
170 225 185 265 
= 


EditT extltem Enabled 
195 85 210 125 
-10 


EditTextltem Enabled 
195 225 210 265 
10 


EditTextltem Enabled 
230 135 245 165 
50 


StatText Disabled 
20 10 40 255 
One Parameter Root Locus Plot Data 


StatText Disabled 
65 25 80 100 
AMin Gain 
StatText Disabled 
65 180 80 260 
AMax Gain 


StatText Disabled 
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170 25 185 80 
X Min 


StatText Disabled 
170 165 185 220 
XMax 


StatText Disabled 
195 25 210 80 
Y Min 


StatText Disabled 
195 165 210 220 
YMax 


StatText Disabled 
230 25 245 125 
Points To Plot 


fee DLOG 


,258 (36) 
Two Parameter Root Locus Plot Data 
30 86 330 426 
Visible NoGoAway 
, 
0 
258 


TYPE DITL 
,208 (36) 
35 


Btnltem Enabled 
35 270 55 330 
Cancel 


Btnitem Enabled 

10 270 30 330 

Plot 

Radioltem Enabled 
120 40 135 265 
Linear Point Interval 


Radioltem Enabled 
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140 40 155 265 
Logarithmic Point Interval 


Radioltem Enabled 
250 130 265 180 
Start 


Radioltem Enabled 
250 185 265 230 
End 


Radioliem Enabled 
250 235 265 288 
Right 


Radioltem Enabled 
250 293 265 338 
Left 


Radioltem Enabled 
270 130 285 180 
Start 


Radioltem Enabled 
270 185 285 230 
End 
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